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ABSTRACT 

We have measured the luminous AGN population in a large sample of clusters of galaxies and find evidence 
for a substantial increase in the cluster AGN population from z ~ 0.05 to z ~ 1.3. The present sample now 
includes 32 clusters of galaxies, including 15 clusters above z = 0.4, which corresponds to a three-fold increase 
compared to our previous work at high redshift. At z < 0.4 we have obtained new observations of AGN 
candidates in six additional clusters and found no new luminous AGN in cluster members. Our total sample of 
17 low-redshift clusters contains only two luminous AGN, while at high redshifts there are 18 such AGN, or 
an average of more than one per cluster. We have characterized the evolution of luminous X-ray AGN as the 
fraction of galaxies with Mr < Mg(z)+ 1 that host AGN with rest-frame, hard X-ray [2-10 keV] luminosities 
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at a median z = 0.72. Our best estimate of the evolution is a factor of eight increase to z = 1 and the statistical 
significance of the increase is 3.8<r. This dramatic evolution is qualitatively similar to the evolution of the star- 
forming galaxy population in clusters known as the Butcher-Oemler effect. We discuss the implications of this 
result for the coevolution of black holes and galaxies in clusters, the evolution of AGN feedback, searches for 
clusters with the Sunyaev-Zel'dovich effect, and the possible detection of environment-dependent downsizing. 
Subject headings: galaxies: active - galaxies: clusters: general - galaxies: evolution - X-rays: galaxies - 
X-rays: galaxies: clusters - X-rays: general 
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1. introduction 

The demographics of AGN in clusters of galaxies have im- 
portant implications for the growth of the supermassive black 
holes at the centers of cluster galaxies, the nature of AGN 
fueling, and the impact of AGN on the intracluster medium 
(ICM) over cosmic time. The luminous, massive elliptical 
galaxies that dominate the galaxy population in the richest 
clusters are also expected (and in some ca ses are measured: 
iHoughton et al.ll2006t iGebhardt et~ai1l2007l) to have the most 
massive black holes in the local universe. As the stars in 
these galaxies appear to have an earlier mean formation epoch 
than those in field galaxies (e.g. Ivan Dokkum & Franxlll996t 
iKelson et~aT1ll997|), the apparent coe volution of black holes 
and galaxies (e.g. Hopkins et al. 2006, and references therein) 
implies that the bulk of their present black hole mass was also 
accreted at earlier times. 

This scenario is also motivated by observations of local 
clusters that clearly show their galaxy populations are more 
qu iescent than local field galaxies. An early demonstration 
by lOsterbrockl dl 960) showed that cluster ellipticals were far 
less likely to have [O II] A3727 emission than field ellipti- 
cals, a result that has since been confirmed by many studies 
(e.g. lGisledfT978l:lDressler et al.lll985l Il999h . One big ques- 
tion that has motivated this work is: Why are galaxy popu- 
lations different in clusters? Numerous physical mechanisms 
have been invoked to explain the relative lack of star forma- 
tion in cluster galaxies, a s well as their higher fraction of el- 
liptical and SO galaxies (lDresslerl[l 980) and relative lack of 
cold gas (e.g. lGiovanelli & Havnesl l985). These include ram- 
Current Address: Department of Astronomy, University of Virginia, P.O. 
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pressure stripping by the ICM dGunn & Gottlll972h . evapora- 
tion of a galaxy's interste llar medium (ISM) by the hot ICM 
(Cowie & Song ailal 19771). tidal effects with the cluster poten- 
tial dFarouki & Shapirolll98lHMerrittlll983HByrd & Valtonenl 
1990), the absence of newly-accreted cold gas (lLarson et al.1 
1980), an d galaxy harassment and mergers dRichstonell 19761 
Moore et all 19961) . 

All of these physical effects may also be important for fuel- 
ing accretion onto the central black holes in galaxies because 
they impact either the available gas supply in a galaxy, an- 
gular momentum transport, or both. The best and perhaps 
only candidate process for fueling the m ost luminous AGN is 
the m erger of two gas-rich galaxies (e.g. lBarnes & Her nquist 
1 19921) and the relative lack of both cold gas and major merg- 
ers is a reasonable explanation for the nearly complete ab- 
sence of QSOs hosted by cluster galaxies. For less luminous 
AGN the case is less clear because an increasing number of 
physical processes such as minor mergers, galaxy harassment, 
various typ es of bars, s tellar mass loss, etc. could also play a 
role (see Martini 200J, for a review). If mechanisms such as 
galaxy harassment and stellar mass loss are important for fu- 
eling low-luminosity AGN, then comparable numbers of low- 
luminosity AGN may be present in clusters and the field. 

Recent studies of the AGN fraction as a function of en- 
vironment with emission-line galaxies from the Sloan Digi- 
tal Sky Survey (SDSS) find that the most luminous AGN are 
rarer in denser environm ents (SDSS; Kauff mann et al.ll2004t 
Popesso & Biviano 2006), although these studies do not sam- 
ple the densest regions of clusters well. This d ecrease is in 
contrast to both lower-luminosi ty AG N in S PSS dMiller et all 
120031) and radio observations dBestl l2004t iBest et al.1 120051) . 
which show that the radio AGN fraction does not decrease 
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significantly in denser environments. X-ray observations with 
Chandra show that the X-ray AGN fraction is larger than ex- 
pected from AGN selection via visible-wavelength emission- 
lines. In previous work we showed that X-ray observations 
identified approximatel y five times as many AGN a s selection 
at visible-wavelengths dMartini et alJl2002l 12006). although 
the precise value of the X-ray excess depends significantly 
on the relative sensitivity and luminosity threshold of the ob- 
servations. This spectroscopic study of X-ray counterparts 
confirmed the many previous studies that suggested a higher 
X-ray AGN populat ion in clusters from surface density ar- 
guments alone (e.g. ICappi et alJl200U ISun & Murray! 120021: 
Ruderman & Ebeling 2005), yet it is still not clear if the X- 
ray AGN fraction is higher than the field value. To date there 
is only weak evidence that the X-ray AGN fraction in clus- 
ters is comparable to the fraction in field early-type galaxie s 
dLehmer et al]l2007llslvakoff et al.ll2008tlArnold et alJl2009h . 
One of the virtues of the emission-line galaxy studies as a 
function of environment is that they can directly calculate the 
fraction of a given galaxy population that hosts AGN as a 
function of environment, even though this technique appears 
to systematically miss AGN in the densest regions relative to 
X-ray and radio selection. 

In addition to a local comparison between AGN in different 
environments, measurement of the evolution of the AGN pop- 
ulation in clusters can constrain the formation epoch for their 
supermassive black holes and the extent of their coevolution 
with the cluster galaxy population. The key early work on 
the evolutio n of galaxies in clusters was by Butcher & Oemler 
dl978l[l984h . who observed a substantial increase in the frac- 
tion of blue galaxies in higher-redshift clusters. The Butcher- 
Oemler effect is interpreted as an increase in the amount of 
star formation and has been confirmed by many other in- 
dicators, in particula r [O II] emission-line galaxy fractions 
(Poggi anti et alj20"06l) and an increase in the numbe r of 24 /xm 
sources in Spitzer obs ervations of distant clusters dBai et al.l 
120071: Saintong e et al J 120081) . The observed increase brings 
the star formation rate (SFR) in cluster galaxies closer to 
those in the field. At a redshift of z ~ 1 and higher, obser- 
vations with Spitzer even find that galaxies in denser environ- 
ments have higher star f ormation rates than lower-density re- 
gions (Elb az et al.l2007l) . which is opposite the trend observed 
in the local univers e. Similar results have also been found 
with deep UV data dHeinis et a l. 2007). The situation is less 
clear when star for mation is measured wi th the [O II] emission 
line because while lPoggianti et all (12008) find that star forma- 
tion do es not strongly depend on environment, ICooper et alj 
(2008) find the specific star formation rate has a similar de- 
pendence on environment at z = and z = 1 , although the total 
star formation rate is higher in clusters at z = 1 than in the 
field. 

The existence of the Butcher-Oemler effect and the many 
indirect arguments outlined above for a connection between 
star formation and black hole accretion suggest that there 
should be an increased AGN population in high-redshift clus- 
ters. An earl y study of the h i gh-red shift cluster 3C295 
at z = 0.46 by iDressler & Gunnl £1983) found evidence for 
three AGN and was an indication that this may be the case; 
however, their relative scarcity precluded a detailed statisti- 
cal study or targeted studies to deliberately identify cluster 
AGN. This situation changed dramatically with the launch 
of Chandra, whose superb sensitivity and angular resolu- 
tion produced a dramatic increase in efficiency for searches 
for AGN, particularly lower-luminosity sources. Just as the 



case for local clusters, Chandra observations of distant clus- 
ters have revealed substantial populations of point sources 
dCappelluti et al.ll2005t iGilmour et al.ll2009l) . Spectroscopic 
confirmation that these point sources are associated with clus- 
ter members has bee n more c hallenging ( Johnson et al.ll2003l : 
iDemarco et al.1 120051) . but in Eastman et al. (2007|) we com- 
bined new observations of MS2053.7-0449 (z = 0.58) with 
archival data on three additional, z > 0.5 clusters and found 
an approximately order of magnitude increase in the fraction 
of Mr < -20 mag galaxies that hosted AGN more luminous 
than L x ,h > 10 43 erg s" 1 in the hard X-ray band (2-10 keV) 
relative to the samp le of ten low-redshift z < 0.32 clusters in 
iMartini et al.l ([2007). These results have since been strength- 
ened with detailed studies of clusters at z ~ 1 with XMM 
dvan Breukelen et al.l2009t) and me asurements of surface den- 
sity excesses in clusters to z ~ 1.5 (Gal ametz et al.ll2009l) . 

In addition to their application to the coevolution of black 
holes and galaxies, an increase in the AGN fraction in clus- 
ters may also impact the ICM. At low redshifts many studies 
have shown that AGN feedback is a viable explanation for the 
absence of substantial reservoir s of cold gas at the centers o f 
clusters (for a recent review see lMcNamara & Nulsenll2007h . 
This feedback is ascribed to AGN associated with the central 
cluster galaxy, which is almost invariably a luminous radio 
source. In our studies of X-ray AGN this is almost the only 
cluster galaxy in which we are insensitive to the presence of 
an AGN because it is challenging to measure even a bright nu- 
clear point source when juxtaposed with the extended emis- 
sion from the ICM that often peaks near the central cluster 
galaxy. Nevertheless, the evolution of AGN in other cluster 
galaxies is likely to be connected to the evolution of the cen- 
tral AGN as the stars in the most luminous cluster galaxies 
have comparable ages. An increase in the net energy pro- 
duction by AGN in higher-redshift clusters is of interest be- 
cause energy input during cluster formation has been invoked 
as an explanation for the minimum e ntropy level in the ICM 
dKaiserj|199lt lEvrard & Henrvlll991l) . AGN remain perhaps 
the most via ble mechanism, if only because most others ca n 
be ruled out (iKravtsov et al.ll2004tfConrov & Ostrikeril2008h . 
although the details of how AGN feedback couples to the 
ICM remain uncertain. Outside of the central galaxy, an in- 
crease in the number of other AGN associated with clusters of 
gala xies may also affect measurement of other clu ster proper- 
ties dBranchesi et al.l 2007; Bignam ini et al.l20 08). Finally, an 
analogous increase in the radio-loud AGN population in high- 
redshift clusters may contamin ate searches for clusters via 
the Sunyaev-Zel'dovich effect (Sunvaev & Zel dovichlll970T) 
at mm and cm wavelengths. As many searches for clusters 
that exploit this effect are in progress, it is important to char- 
acterize the potential impact of evolu tion of the cluster A GN 
population on these experiments (e.g. lLin & Mohr 20Q3). 

In the next section we describe our expanded high-redshift 
data, as well as the selection criteria for X-ray AGN we em- 
ploy at all redshifts. We then describe our new observations of 
low-redshift clusters in ^3] These two datasets are combined 
to calculate the cluster AGN fraction and its evolution in $4] 
followed by an examination of the properties of the cluster 
AGN in |5] We discuss the implications of these results, par- 
ticularly on the coevolution of black holes and galaxies, in 
$6] and conclude with a summary of our results. Throughout 
this paper we assume that the cosmological parameters are: 
(n M ,n A ,h) = (0.3, 0.7, 0.7) where H = lOO/i km s^Mpc" 1 . 
All absolute magnitudes quoted in this paper assume h = 0.7. 
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2. HIGH-REDSHIFT DATA 

Two large surveys have obtained redshifts for substantial 
numbers of galaxies with X-ray counterparts in many deep, 
archival Chandra observations that include substantial num- 
bers of high-redshift clusters of galaxies. These are the 
Serendipitous E xtragalactic X-ray Source Identification Pro- 
gram (SEXSI: lHarrison et alj|2003t lEckart et al. 2005, HH) 
and the Chandra Multiwavelengfh Project (ChaMP ; lKim et alJ 
I2004allbl: iGreen et alj|2004t ISilverman et~aT]|2005l) . We have 
investigated the fields surveyed by both SEXSI and ChaMP to 
identify datasets that contain clusters of galaxies with z > 0.4 
and have sufficient depth to identify Lx.h > 10 43 erg s" 1 (rest 
frame 2-10 keV) AGN at the cluster redshift. 

The SEXSI survey published sp ectroscopic redshif ts for 
27 archival Chandra observations in lEckart et alJ (120061) that 
were selected to identify hard X-ray sources over the flux 
range of 10~ 13 - 10~ 15 erg s _1 cm~ 2 and isolate those respon- 
sible for the hard X-ray background. The specific selection 
criteria for the fields were that they must be high Galactic lat- 
itude (\b\ > 20°) and be obtained with either the I or S modes 
of the Advanced Camera for Imaging Spectroscopy (ACIS; 
iBautz et alii 19981) when no grating was used. The X-ray lu- 
minosities quoted by SEXSI are based on spectral fits that 
assume a V = 1 .5 power law and intrinsic absorption Nh at the 
source redshift, although they quote the observed luminosi- 
ties (not corrected for obscuration) and provide the best-fit 
Nh value. The average spectroscopic completeness is 67% 
(see {34.2l below) for sources with R < 24.4 mag on the Vega 
system. Nine of the 27 SEXSI fields include clusters of galax- 
ies with z > 0.4 and we include seven 2 in our sample. As one 
field contains 3 clusters, we list nine clusters from SEXSI in 
Table [TJ 

The ChaMP survey published spe ctroscopic redsh i fts for 
20 archival Chandra observations in ISilverman et al.l (2005) 
that were similarly selected for depth, high Galactic latitude 
(\b\ > 20°), and no special observing modes. The spectro- 
scopic completeness of ChaMP is 77% at r ' < 22.5 mag, 
where r' is on the SDSS photometric system dFukugita et al.l 
1996, and r' AB = R Vega + 0.11). Their X -ray luminosities are 
based on spectral fits that assume a T = 1 .9 power law and 
intrinsic absorption Nh at the source redshift, as well as the 
appropriate Galactic absorption, although they also quote the 
observed luminosities (only correct ed for Galactic absorp- 
tion). The final sample presented in ISilverman et al.1 (2005) 
was restricted to X-ray sources with L x > 10 42 keV in the 2-8 
keV band in order to insure all are AGN. Most (69%) are spec- 
troscopically classified as broad-line AGN (BLAGN). Nine 
of these 20 ChaMP fields include clusters of galaxies with 
z > 0.4 and we include eight 3 of these in our study (see Ta- 
bleQ}. Two of these clusters are common to both ChaMP and 
SEXSI (MS2053.7-0449 and RXJ1716.4+6708) and therefore 
the final sample has fifteen clusters with z > 0.4. While spec- 
troscopic dat a for X-ray sources in other high-redshift clus- 
ters exist (e.g. Johnson et al. 2006), we limit our high-redshift 
sample to these fifteen to maximize the uniformity of the 
dataset. 

2 RX J1350. 0+6007 was not targeted for spectroscopy and the X-ray data 
for CL0442+0202 (z = 1.11) were sufficiently shallow (/ = 44ks) that they 
may not be complete to Lx.h = 10 43 erg s _1 . In addition, Stern et al. ( 2003) 
classify CL0442+0202 as an overdensity that has not yet collapsed, rather 
than as a cluster. 

3 We exclude CL J0152. 7-1357 (z = 0.831) because the exposure time is 
shorter than the others at / = 34.6 ks and therefore the X-ray data may not be 
complete to Lx.h = 10 43 erg s _1 . 



We have also compiled additional data for each cluster 
listed in Table [TJ that will be important for our subsequent 
analysis. One quantity is the center of the cluster, which 
is needed to determine if a given AGN falls within the pro- 
jected virial radius of the cluster. We associate the center of 
each cluster with the centroid of the extended X-ray emis- 
sion. While these coordinates do not always agree with the 
standard coordinates quoted in the literature, this assumption 
makes our analysis more uniform. The redshift and velocity 
dispersion are also needed to determine if an AGN is within 
the cluster. In most cases velocity dispersions for these clus- 
ters are available in the literature and we quote the origin 
of the measurement we adopt in the table. When the veloc- 
ity dispersion has not been measured, we estimate this quan- 
tity fr om the X-ray t emperature and the a-Tx relationship 
from IXue~& Wu (2000). Specifically, we used the relation 
cr = io 2 -5i±ooi r o.6i±o.oi ^ s -i derived from meir combined 

group and cluster sample with orthogonal distance regression 
(Feigelson & Babul ll992l) . Based on their data, we estimate 
that there is a 30% uncertainty in a at fixed T. 

One potentia l conc ern for our subsequent analysis is that 
the IXue & Wul (120001) a-T relation may not hold at higher 
redshift. iLubin et al.1 (120041) investigated this point for several 
optically-selected clusters and found that they were 2-9 times 
cooler than expected from the local relation; however, the dif- 
ference was much less stark for X-ray selected, high-redshift 
clusters simil ar (and in several c ases identical to) those pre- 
sented here. Fang et al. (2007) showed that high-redshift, 
X-ray selected clusters are consistent with the low-redshift 
Lx - a relation, although spectroscopically-s elected groups 
and c lusters do not agree as well (see also lAndreon et al.1 
l2008h . 

Finally, we have calculated the projected s i ze of the 
virial radius for each cluster following iTreu et ail (120031) and 
throughout this paper we associate the virial radius with R200, 
the radius within which the cluster is a factor of 200 over- 
density. Of the thre e clusters we have in common with 
iPoggianti et all d2006l) . for 3C 295 and MS 1054-03 we adopt 
nearly the same a and our R 200 estimate is nearly identical to 
theirs, while for MS0015. 9+1609 we adopt a slightly larger 
velocity dispersion (1234 km s -1 from Carlberg et al.l (1 19961) 
rather than their 984 km s" 1 ) and consequently infer a larger 
radius. 

Because the most relevant ChaMP measurements are the 2- 
8 keV luminosity, rather than 2-10 keV luminosity, we multi- 
ply the ChaMP 2-8 keV luminosities by a factor of 1.2. This 
correction factor was calculated for a T = 1 .7 power law with 
PIMMS. There is some uncertainty in this correction factor 
because not all AGN have this power-law form, particularly 
as we assume this correction for their observed rather than in- 
trinsic (unobscured) spectra, but this is not a significant effect 
compared to other sources of systematic errors that we discuss 
below. There are no additional AGN from ChaMP that enter 
the sample after this step because there are none just below the 
10 43 erg s -1 threshold in the 2-8 keV band. We also estimated 
the difference in luminosity for an AGN calculated with the 
T = 1 .5 power law employed by SEXSI, the T = 1.9 employed 
by ChaMP, and a T = 1 .7 power law to determine if these dif- 
ferences would cause any sources to fall in or out of the same 
and none would do so. In the two clusters observed by both 
ChaMP and SEXSI, there is one cluster AGN common to both 
surveys: CXOSEXSI J171636.9+670829. The redshifts from 
the two surveys agree exactly (z = 0.795) and the luminosities 
agree well: Lz.2-10 = 10 44 erg s"'and Lx.i-ss = 10 43,88 erg s" 1 . 
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TABLE 1 
High-Redshift Cluster Sample 



Cluster 


etc 


5c 


z 


(j [km/s] 


o-Ref 


T x [keV] 


7>Ref 


R 2 oo [Mpc] 


Spectra 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


MS 1621. 5+2640 


16:23:34.9 


+26:34:21 


0.43 


735 


1 


7.6 


1 


1.42 


SESXI 


3C 295 


14:11:20.5 


+52:12:09 


0.46 


1642 


1 


5.3 


1 


3.12 


SESXI 


MS 045 1. 6-0305 


4:54:11.1 


-03:00:55 


0.538 


1371 


2 


8.1 


1 


2.49 


ChaMP 


Mo 0UlJ.y+16UV 




+ 10.ZO.UO 


U.j4I 




z 




z 


Z.Z4 


cnaMr 


RX J0848.7+4456 


8:48:47.6 


+44:56:16 


0.574 


670 


3 


3.2 


3 


1.19 


SESXI 


MS 2053.7-0449 


20:56:21.3 


-04:37:49 


0.583 


865 


4 


5.2 


1 


1.53 


SESXI.ChaMP 


RX J0542.8-4100 


5:42:49.8 


-41:00:07 


0.634 


1101 


3 


7.9 


3 


1.89 


ChaMP 


RX J2302.8+0844 


23:02:48.3 


+08:43:48 


0.722 


993 


3 


6.6 


3 


1.61 


ChaMP 


MS 1137.5+6625 


11:40:22.1 


+66:08:14 


0.782 


967 


3 


6.3 


1 


1.52 


ChaMP 


RXJ1317.4+2911 


13:17:22.0 


+29:11:24 


0.805 


531 


3 


2.2 


1 


0.82 


SESXI 


RX J17 16.4+6708 


17:16:49.3 


+67:08:25 


0.813 


1445 


1 


6.6 


1 


2.22 


SESXI.ChaMP 


MS 1054-03 


10:56:55.7 


-03:37:39 


0.831 


1156 


5 


7.8 


1 


1.76 


ChaMP 


RDCS J09 10+5422 


9:10:44.7 


+54:22:04 


1.11 


675 


6 


3.5 


1 


0.87 


SESXI 


Lynx E 


8:48:58.3 


+44:51:51 


1.261 


740 


7 


3.8 


4 


0.88 


SESXI 


Lynx W 


8:48:34.2 


+44:53:35 


1.27 


650 


8 


1.7 


4 


0.77 


SESXI 



NOTE. — Cluster sample and properties derived from the present study. Columns are: (1) Cluster name; (2 and 3) RA and DEC for the centroid of the extended X-ray emission; 
(4) redshift; (5) velocity dispe rsion; (6) reference for the velocity dispersion; (7) X-ray temperature in keV ; (8) reference for th e X-ray temperature; (9) estimate of the virial radius 
in Mpc (e.g., Treu et al. 200 J); (10) origin of m ost of the spect r a. Ref ere nces for velocity d ispe rsion are: 1 : IGir ardi & Mezzetti 1 2001); 2: Carlberg et al. (1996); 3: d erived from the 
X-ray temperature following lXue & Wul 120001) ; 4: iTran et all 120051); 5 : [ftan et al. ] 2007); 6: Mei et al. r200fl) ; 7: fro m weak lensing es timate lJee et alj 120061) ; 8: IStanford et"al] 
t200ll) . References for X-ray temperatures are: 1 : 1 Vikhlinin et al J 120021) ; 2: lEbelins et al. (2007); 3: Ettori et al. 1 2004); 4 : |jeeetal] 120061) . 



TABLE 2 

High-Redshift Cluster AGN Sample 



AGN 


Cluster 


z 


R [mag] 


logLx.H [ergs '] 


Sv/a 


AS [arcmin] 


R/R200 


Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


CXOSEXSI J141 127.4+521 131 


3C295 


0.451 


19.78 


43.4 


1.13 


1.23 


0.14 


ALG 


CXOSEXSI J141 123.4+521331 


3C295 


0.472 


19.05 


43.8 


1.5 


1.45 


0.16 


BLAGN 


E0015+162 


MS0015.9+1609 


0.553 


18.41 


45.48 


1.89 


3.35 


0.58 


BLAGN 


CXOSEXSI J084858.0+445434 


RX J0848.7+4456 


0.573 


19.58 


43.8 


0.28 


2.5 


0.83 


BLAGN 


CXOMP J054248.2-410140 


RDCSJ0542-4100 


0.634 


20.64 


43.24 





1.58 


0.32 


NELG 


CXOMP J05425 1.4-41 0205 


RDCSJ0542-4100 


0.637 


19.63 


43.35 


0.5 


1.99 


0.33 


ALG 


CXOMP J054259.5-410241 


RDCSJ0542-4100 


0.638 


20.50 


43.37 


0.67 


3.16 


0.63 


NELG 


CXOMP J054240.8-405626 


RDCSJ0542-4100 


0.639 


20.89 


43.67 


0.83 


4.05 


0.81 


NELG 


CXOMP J054255.0-405922 


RDCSJ0542-4100 


0.644 


22.08 


43.08 


1.67 


1.24 


0.25 


NELG 


CXOMP Jl 14022.0+660816 


MS 11 37+6625 


0.786 


20.37 


43.24 


0.7 


0.04 


0.01 


BLAGN 


CXOSEXSI J17 1636.9+670829 


RXJ17 16.4+6708 


0.795 


22 


44 


2.06 


1.19 


0.24 


ELG 


CXOSEXSI J131718.8+291 111 


RX J1317.4+291 1 


0.803 


21.98 


43.3 


0.63 


0.68 


0.38 


BLAGN 


CXOSEXSI J171703.8+670900 


RXJ17 16.4+6708 


0.812 


21.79 


43 


0.11 


1.53 


0.31 


ELG 


CXOSEXSI J171714.5+671 136 


RXJ17 16.4+6708 


0.815 


22.68 


43.2 


0.23 


4.02 


0.82 


ELG 


CXOMP J 105650.6-033508 


MS 1054-03 


0.818 


21.76 


43.22 


1.84 


2.82 


0.73 


BLAGN 


CXOU J091043. 3+542152 


RDCSJ09 10+5422 


1.104 


24 


43.06 


1.26 


0.29 


0.16 


AGN2 


CXOSEXSI J084905.3+445203 


LynxE 


1.266 


24.61 


43.8 


1.11 


1.27 


0.74 


ELG 


CXOSEXSI J08483 1 .6+445442 


LynxW 


1.267 


25.42 


43.2 


0.61 


1.23 


0.8 


ELG 



NOTE. — AGN in high-redshift clusters of galaxies. Columns are: (1) AGN name; (2) Cluster; (3) AGN redshift; (4) fi-band magnitude; (5) Rest-frame, hard-X-ray luminosity 
(2-10 keV); (6) Velocity offset from the cluster systemic velocity normalized by the cluster velocity dispersion; (7) Projected radial offset relative to the centr oid of the X-ray gas i n 
arcminutes; (8) Projected radial offset normalized by the cluster virial radius; (9) Spectroscopic classification. The tf-band magnitude of E0015+162 is from lOrndahl et alj (2003). 
The remaining values are from either Eckart et al. 1 2006) for the SEXSI sample or from Silverman et al. 1 2005) for the ChaMP sample (although corrected from r' to R as noted in 
Section|2). The 2-8 keV X -ray luminosities from Silverman et al. 1 2005) have been corrected to the 2-10 keV band as described in Sectionl2l 



We also correct the ChaMP r' measurements to the Vega R 
band as discussed above. Based on the magnitudes of these 
sources and a simple ^-correction, we estimate that none of 
these sources falls below our galaxy luminosity threshold. As 
these are fairly luminous AGN, in some cases the AGN may 
dominate the total flux and we may have overestimated the 
host galaxy luminosity. E0015+162 dMargon et al.ll 1983D is 
the most X-ray luminous AGN in our sample by over an order 
of magnitude and is a useful case study to test the importance 
of this concern. This AGN has a total R = 18.41 mag and a 
fainte r host galaxy magnitude of R = 19.8 mag (Orndahl et al. 
2003), which corresponds to a factor of 3.6 in flux. If the 
other AGN have similar or smaller Lr/Lx ratios (such as due 
to obscuration), then we expect their AGN contribution to the 
measured /?-band flux to be negligible because they are all 
much less luminous than E0015+162. 



We identify AGN in these clusters with the following four 
criteria: 1) The hard X-ray luminosity must be Lx,h > 10 43 
erg s _1 ; 2) The AGN redshift must fall within 3<r of the clus- 
ter mean redshift, where a is the cluster velocity dispersion; 
3) The AGN must fall within the projected virial radius /?200 
of the cluster; 4) The absolute magnitude of the host galaxy 
must be greater than Mr =M%(z) + \ mag. M ost of these cri- 
teria were adopted from Eastman et al. (2007), although the 
absolute magnitude criterion i s dif ferent and we discuss our 
motivation for this choice in ^43] below. With these criteria 
we identify 18 AGN in the 15 clusters with z > 0.4, or an aver- 
age of more than one per cluster. The properties of the z > 0.4 
AGN are presented in Table [2] 



3. NEW LOW-REDSHIFT OBSERVATIONS 
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TABLE 3 
New Low-Redshift Clusters 



Cluster 


Ct c 


So 


z 


cr [km/s] 


cr Ref 


T x [keV] 


T x Ref 


i?,oo [Mpc] 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 



Abell 1240 


11:23:37.3 


+43:06:54 


0.1590 


698 


1 






1.64 


Abell 1942 


14:38:22.0 


+03:40:07 


0.2240 


903 


2 


5.6 


1 


1.96 


Abell 2125 


15:41:13.2 


+66:16:01 


0.2465 


1113 


3 


3.2 


2 


2.39 


MS1455.0+2232 


14:57:15.1 


+22:20:29 


0.2578 


1032 


4 


5.5 


3 


2.20 


ZwCl 1358.1+6245 


13:59:50.6 


+62:31:04 


0.3280 


1003 


4 


6.5 


3 


2.06 


MS 15 12.4+3647 


15:14:22.4 


+36:36:21 


0.3720 


575 


4 


3.6 


3 


1.15 



NOTE. — New low-redshift clusters and their properties derived from the present study. Columns are: (1) Cluster name; (2 and 3) RA and DEC for the centroid of the extended 
X-ray emission; (4) redshift; (5) velocity dispersion; (6) reference for the velocity dispersion; (7) X-ray temperature in keV; (8) ref erence for the X-ra y temperature; (9) estimate 
of the virial radius in Mpc iTreu et al. 2003). References for velocity dispersion are: 1: derived from the X-ray luminosity following Xue & Wu ( 2000); 2: derived from the X- ray 
temperature followinglXue & Wul <200fJ) : 3: IMiller et all (20041) : 4: IBorgani et alU1999l) . References for X-ray temperatures are: 1 : IQta & Mitsud A 12004) ; 2: IWang et aT] 12004) ; 3: 
iMushotz kv & Scharf 1 19971) . 



TABLE 4 
Chandra OBSERVATION LOGS 



Cluster 


OBSID 


Detector 


T 


Lx,H.Um 








(ks) 


(10 41 ergs _1 ) 


(1) 


(2) 


(3) 


(4) 


(5) 


Abell 1240 


4961 


ACIS-I 


51.3 


1.2 


Abell 1942 


3290 


ACIS-I 


57.5 


2.2 


Abell 2125 


2207 


ACIS-I 


81.5 


1.9 


MS 1455.0+2232 


4192 


ACIS-I 


91.9 


1.8 


ZwCl 1358.1+6245 


516 


ACIS-S3 


53.0 


2.6 


MS 1512.4+3647 


800 


ACIS-S3 


36.4 


4.6 


NOTE. — Chandra Observation Log. Columns 


are: (1) Cluster 



targeted; (2) Observation ID of Chandra data; (3) Detector used; (4) 
Usable exposure; (5) Estimate of the 2.0-8.0keV luminosity limit 
of the observation for a cluster galaxy. 

AGN more luminous than Lx.h = 10 43 erg s" 1 are suffi- 
ciently rare in low-redshift clusters that Poisson uncertainties 
(as opposed to sources of systematic errors) from the low- 
redshift sample may dominate the statistical significance of 
any evidence of evolution. Our previous study of ten clus- 
ters with z < 0.32 only identified one AGN above this lu- 
minosity threshold (Mar tini et al.ll2006l) . while our more re- 
cent observations of three additional clusters (all at z < 0.08) 
hav e identified only one additional AGN above this luminos- 
ity dsTyakoff et al. 2008). We have therefore studied six addi- 
tional clusters with 0.15<z<0.4to find other X-ray AGN 
more luminous than Lx,h = 10 43 erg s" 1 with a combination 
of Chandra archival data and follow-up spectroscopy of can- 
didate cluster X-ray AGN at the MDM Observatory. These 
clusters were selected to be the nearest massive clusters in the 
Chandra archive whose estimated virial radii fit within the 
Chandra ACIS field of view (FOV) and were accessible dur- 
ing our observing runs. The new clusters and their physical 
properties are listed in Table [3] 



3.1. Chandra X-ray Analysis 

The X-ray observation s were processed follo wing the same 
techniques employed by Siva koff et al.l d2008l) . We reduced 
all data using CIAO 3.4 4 with CALDB 3.3.0.1 and NASA's 
FTOOLS 6.0 5 . The observations are summarized in Table @] 
Only minor differences in reduction were required for these 
archival observations. The majority of the clusters had data 
with an aimpoint centered on the four ACIS-I chips (~ 17' 
FOV) and frame times of 3.1 s. These data were telemetered 
and cleaned in Very Faint mode. The more distant clusters, 
ZwCl 1358.1+6245 and MS 1512.4+3647, were observed 



with the aimpoint placed on the ACIS-S3 detector (8.4' FOV) 
and had frame times of 3.3s. Their data were telemetered 
and cleaned in Faint mode, and thus have a slightly higher 
background. As all observations were operated at -120° C 
the X-ray data were corrected for the time dependence of 
the gain and the charge-transfer inefficiency with their pho- 
ton energies determined using the gain file acisD2000-01- 
29gain_ctiN0006.fits. The archival data of all observations 
already had applied the newest tools to detect hot pixels and 
cosmic ray afterglows. We only consider events with ASCA 
grades of 0, 2, 3, 4, and 6. Known aspect offsets were applied 
for each observation. All observations were corrected for 
quantum efficiency degradation and had exposure maps deter- 
mined at 1 .5 keV. We excluded bad pixels, bad columns, and 
columns adjacent to bad columns or chip node boundaries. 
We also filtered out times when the blank-sky rate was more 
than three times the expected blank-sky rate derived from cal- 
ibrated blank-sky backgrounds to avoid the most extreme pe- 
riods of high background ("background flares") that Chandra 
may encounter. MS 1512.4+3647 had two separate pointings 
and this introduced difficulties into our standard processing. 
We therefore excluded the shorter second pointing, which ac- 
counted for less than 25% of the total integration time. 

To detect X-ray sources that are potential X-ray AGN in 
these clusters, we applied the wavelet detection algorithm 
(CIAO WAVDETECT) with scales ranging from 1 to 64 pixels 
in steps of \/2 factors and required a source detection thresh- 
old of 10~ 6 . Source detection was only performed in regions 
with an exposure of greater than 10% of the total for the ob- 
servation. Our source detection threshold corresponds to < 4 
falsely detected X-ray source s (due to a statist ical fluctuation) 
for each observation. Using iKim et ail (120071) . we have es- 
timated the statistical X-ray positional uncertainty (la) due 
to WAVDETECT. In Tabled we list an estimated limiting X- 
ray luminosity for each observation that corresponds t o five 
counts on axis (for consistency with lMartini et al.ll2.Q06t) . For 
our analysis we concentrated on sources with at least 20 broad 
(0.3-8.0 keV) X-ray counts. These sources are unlikely to be 
due to statistical fluctuations except where they are coincident 
with ICM emission. 

We used ACIS Extract 3.131 6 to create source extraction 
regions enclosing 90% of the flux in the X-ray PSF and to de- 
termine a masking radius that encircled 97% of the flux. For 
most of the sources, whose photons had median energies of 
~ 0.6-2.6keV, we determined the regions assuming the PSF 
at 1.497keV. A few sources had harder emission and their 
PSF was calculated assuming an energy of 4.51keV. In a 



http : //asc . harvard . edu/ciao/ 



http : //heasarc .gsf c .nasa . gov/docs/ software/ lheasoft / http://www.astro.psu.edu/xray/docs/TARA/ae_users_guide.html 
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relatively small number of crowded regions, the PSF fraction 
was reduced to prevent overlapping source extraction regions. 
We also used ACIS Extract to correct the (CIAO WAVDETECT) 
position to the mean position of detected events for sources 
within 5'of the observation aimpoint or to the position that 
best correlated with the PSF for sources beyond 5'of the ob- 
servation aimpoint. These new positions were registered with 
an optical catalog from /?-band images (see below) to correct 
the absolute astrometry and determine the absolute astromet- 
ric precision of each Chandra observation (0.3-0.5"). The 
statistical significance of each detection was added in quadra- 
ture with the absolute astrometric precision to estimate the 
total X-ray positional precision. We measured the counts in 
three energy ranges: the broad (0.3-8 keV), soft (0.3-2 keV), 
and hard (2.0-8.0 keV) bands. The observed fluxes in these 
bands were derived assuming a T = 1 .7 power-law spectrum 
with Galactic absorption. We then calculated the rest-frame 
luminosity in the broad band (0.3-8 keV) and the classic hard 
band (2-10 keV) for all sources with redshifts (see 



3.2. MDM Photometry 

7?-band images of these clusters were obtained at the MDM 
Observatory 2.4m Hiltner telescope with the Echelle CCD 
camera during a run from the night of 28 May 2007 to 3 June 
2007. Because the FOV of the CCD camera (~ 9.5' x 9.5') 
is smaller than the ACIS-I FOV (- 17' x 17'), we imaged a 
2x2 mosaic to cover the Chandra area, with each panel con- 
sisting of 3 x 300 s exposures. All images were trimmed, bias- 
subtracted and flat-fielded with the CCDPROC package within 
IRAF 7 . 

Sources were catalog ed with the SExtractor package 
dBertin & Arnoutsl fT996). Aperture magnitudes from these 
catalogs were calibrated with multiple observations of 
standard star fields from the data compiled by P. B. 
Stetson 8 onto the Vega magnitude system. Only data 
from the last night, which includes each quadrant of 
Abell 1240 and ZwCl 1358.1+6245, the north-east quadrant 
of MS 1512.4+3647, and 1 x 300s exposures of each quad- 
rant of Abell 2125, were taken under photometric conditions. 
Our derived photometric solution for this night was precise 
to 0.03 mag. As all of these clusters except for Abell 2125 
were imaged with SDSS, we cross-correlated aperture mag- 
nitudes from all images on this run with stars in the SDSS 
DR5 catalog. After correcting to R (Vega) 9 , our derived 
photometric solution for 3 June, which includes a color cor- 
rection term, is accurate to 0.01 mag and precise to 0.06 
mag. The poorer precision compared to our photometric so- 
lution appears to be only partially due to the dispersion in 
the Vega correction (overlapping sources between quadrants 
of our own observations indicate typical photometric preci- 
sions of 0.05 -0.08 mag). We therefore adopted the SDSS 
cross-calibration technique to photometrically correct all ob- 
servations on non-photometric nights, except for observations 
of Abell 2125. For Abell 2125, non-photometric observations 
were cross calibrated with the single photometric exposures 
for Abell 2125. As we do not have complete multi-band data, 
we report only the magnitudes assuming no color correction. 

7 IRAF is distributed by the National Optical Astronomy Observatory, 
which is operated by the Association of Universities for Research in Astron- 
omy (AURA) under cooperative agreement with the National Science Foun- 
dation, 

8 \protect ittp://cadcwww.hia.nrc.ca/standards 



The exclusion of the color correction term does not signifi- 
cantly decrease the precision of our photometric solutions. 

We calculated astr ometric so lutions for the images with the 
WCSTools package (Mink 2002), package and then produced 
the final, calibrated mosaics with the SWARP 10 package. A 
final source catalog was extracted with SExtractor and used to 
register the astrometry of the X-ray observations. We consider 
only the SExtractor AUTO magnitudes, which is an automatic 
aperture magnitude designed to give precise estimates of total 
magnitudes for galaxies. As nearby, detected neighbors are 
removed and replaced by mirroring the opposite side of the 
aperture where available, these magnitudes are suitable for 
our relatively crowded fields. All X-ray sources that would 
be more luminous than Lx.h = 10 43 erg s" 1 at the cluster red- 
shift that were also associated with galaxies and that would 
be more luminous than M^(z)+ 1 at the cluster redshift were 
then targeted for the highest-priority spectroscopic observa- 
tions, with the exception of sources heavily contaminated by 
ICM emission. We also identified other candidate cluster X- 
ray AGN, specifically those that would have Lx.h > 10 42 erg 
s" 1 , as lower-priority spectroscopic targets. 

3.3. MDM Spectroscopy 

We obtained low-resolution spectroscopy of these candi- 
dates with the 2.4m Hiltner telescope with the CCDS, a Boiler 
& Chivens spectrograph, during a run from the night of 28 
April 2008 to 3 May 2008. The slit widths were determined 
by the nightly seeing conditions and were either 1.0" or 1.5". 
At least two exposures of every candidate were obtained and 
total exposure times varied from 120 s to 9000 s. Five sets 
of internal and twilight flats were taken over the entire run, 
while comparison lamps were observed before and/or after 
every candidate. 

The files were trimmed and bias-subtracted with the 
ccdproc package within IRAF and bad pixels were deter- 
mined from a ratio of flat-field images and were fixed in every 
image. The individual flat-field images from internal lamps 
revealed a complex wavelength and slit-dependent flat-field, 
most likely due to some reflection. To model this complex re- 
sponse, we first median smoothed the internal flat-fields (over 
11x11 pixels) and then Gaussian-smoothed (a = 1 1 pixels) 
over the dispersion axis. The ratio of the internal flat-field 
to the modeled internal flat-field was adopted as the true in- 
ternal flat-field. An illumination correction was then created 
from the twilight flat-fields and applied to make the final set of 
flat-field corrections to remove fringing in the spectra. After 
each spectrum was proper ly flat-fielded, we re jected cosmic 
rays using LA. Cosmic 11 (Ivan Dok kum 2001). A fourth or- 
der wavelength-solution was calculated for each set of HgNe 
comparison spectra, resulting in a typical RMS of ~ 0.1 A 
pixel" 1 . Thereafter, standard aperture extraction of the spec- 
tra was used to remove the night sky emission and produce 
one-dimensional, logarithmically interpolated spectra with a 
dispersion of ~ 3 A pixel" 1 . The spectra extend from approx- 
imately 3650A to 7250A. We extracted both the signal and 
noise for each final spectrum of a source. 

We adapted the Princeton/MIT SDSS Spectroscopy rou- 
tines 12 to calculate redshifts. This technique cross-correlates 
the spectra in pixel space with template spectra, with each 

10 http://terapix.iap. fr/rubrique.php?id_rubrique=49 

11 \protect "ittp://www.astro. yale.edu/dokkum/lacosmic/ 



protect rttp://www.sdss.org/dr7/algorithms/sdssUBVRITransform.html#LuptoM2^5rDtect "ittp://spectro. princeton.edu/idlspec2d_doc. html 
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TABLE 5 

New, Lower-Luminosity Cluster X-ray AGN 



CXOU ID 


z 


zref 


R 


flflag 


fx.s 


fx.H 


fx.B 


Lx.B 


Lx.H 


Xflag 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


(11) 



J135950.5+623 106.3 
J143821.8+034013.3 
J145714.7+221933.6 
J145715.0+222034.5 
J151422.5+363620.7 
J154101.9+661627.1 
J154101.9+661721.4 
J1541 17.3+661923.6 



0.32717 ±0.00038 
0.22479 

0.24852 ±0.00025 
0.25772 ±0.00015 
0.3718 

0.24564 ±0.00045 

0.2567 

0.2465 



17.80 ±0.05 
16.44 ±0.06 

20.04 ±0.07 
16.82 ±0.07 

18.05 ±0.06 
17.19±0.08 
19.36 ±0.08 

18.81 ±0.08 



9.5+ 11 
2.20+Mff 

2 40+8:61 

20 2+ 4 ' 9 

3 Q8+ - 98 
2 7g+0.62 



3.7$ 
21.6!; 
3 2+ 2 ' 5 

21+ - 76 

7 1 +L7 

-1.4 

1 46+' 2 



20.6+ 2 - 2 
4 - 8 -15 

44 4+W 
8 4+ 1 - 9 

17.0+jj 
4 1 5+ 1J 



70.6+™ 
7.2+1' 
11.2+1 
93+ 

34.3+" 
7 6+' : ' 



46.2+5-0 

4 7+ i;7 
^•'-1.4 

7 3+ 1 - 5 
61+ 12 
24.9!" 

5 5+ L2 
22.4+|] 

5 0+" 



NOTE. — Chandra Ob ser vation Log. Co l umns are: (1) Name of X-ra y source; (2) Redshift (3) References for redshift are: 1: this work; 2: SDSS 
lAdelman-Mc Carthv 2008); 3 Abraham et al. ( 1998); 4: IMiller et alj 120041) ; (4) fi-band magnitude; (5) Flags for photometry are: (0) no flag; (1) may be 
contaminated by nearby neighbors or bad pixels; (2) blended with nearby neighbors; (3) both; (6-8) Soft [0.3-2 keV], Hard [2-8 keV], and Broad [0.3-8 keV] 
band flux in the observed frame in units of 10~ 15 ergs~'cirr 2 . (9-10) Broad [0.3-8 keV] and Hard [2-10 keV] band luminosity in the rest-frame in units of 
10 41 ergs" 1 corrected for Galactic absorption. (11) X-ray flags are: (0) no flag; (1) contaminated by ICM peak. Note that CXOU J145715. 0+222034.5 is the 
BCG and we subtracted a multi-component beta model for the ICM to compute the quoted fluxes and luminosities. 



pixel weighted by the inve rse of its vari a nce, a nd is simi- 
lar to the technique used in Marti ni et al.l (120061) . The tem- 
plate spectra include a set of four eigenspectra for galaxies, 
four eigenspectra for quasars, and forty eigenspectra for stars. 
The five best galaxy redshifts for -0.01 < z < 1.00, five best 
quasar redshifts for 0.0033 < z < 7.00, and forty different stel- 
lar redshifts for -0.004 < z < 0.004 are found and ordered by 
the reduced x 2 of their fit. We adopted the best-fit redshift 
and classification for each source. To ascertain the quality of 
the fit and the errors to the redshift, we resampled each spec- 
tra 100 times randomly according to its noise characteristics 
and reran the cross-correlation routine. Both the dispersion 
in best-fit redshifts and the best-fit spectral type were used to 
qualify the spectral classification quality. If the dispersion in 
redshift was relatively low (a z < 0.01), > 68% of the best- 
fit redshifts were within 3cr z of our adopted redshift, and had 
the same spectral type (i.e., galaxy, quasar, or a similar stellar 
type) we consider this a secure redshift. Typically the maxi- 
mum SNR of these spectra were > 5 pixel" 1 . 

We did not identify any AGN in these clusters with 
Lx.h > 10 43 erg s" 1 , although we did identify several lower- 
luminosity AGN in these clusters. Data for the lower- 
luminosity X-ray sources are provided in Table [5] and include 
several sources with spectroscopic measurements from the lit- 
erature. The spectroscopic observations of Abell 1240 and 
MS 15 12.4+3647 are complete for all candidates that would 
have Lx.h > 10 42 if at the cluster redshift, while the other 
four clusters are not complete to this luminosity limit. We 
have also measured redshifts, 7?-band magnitudes, and X-ray 
fluxes and luminosities for numerous additional sources not 
associated with these clusters and their properties are listed in 
Table [6] As for the high-redshift clusters, several of the low- 
redshift clusters do not have direct velocity dispersion mea- 
surements. For Abell 1942 we estimated this quantity from 
the X-ray temperature. For Abell 1240 IXue & Wul d2000l) 
quote kT = 3.83 ke V from jMushotzkv & Scharflll997l) . but 
in fact the value in Mushotzky & Scharfl (1 19971) appears in- 
stead to be for Abell 1242. As we could not identify an- 
other Tx value in the literature , we used the mea surement of 
L bo i = 2.71 x 10 44 erg s" 1 from lDavid et al.| (|1999|) and the re- 
lation a = 10 2 ' 76 L° 19 
the velocity dispersion 



L£ ly derived by lXue & Wul (|2000l) to estimate 



4. CLUSTER X-RAY AGN FRACTION 



We require two quantities to estimate the AGN fraction in 
these clusters: the number of AGN above our hard X-ray lu- 
minosity threshold hosted by galaxies with Mr < Mr(z) + 1 
and the total number of cluster galaxies above this magni- 
tude threshold. For our low-redshift cluster sample, we have 
complete data to our X-ray threshold and reasonably complete 
data for the other cluster galaxies for about half of the clus- 
ters. For the high-redshift sample we have incomplete knowl- 
edge of both quantities. The AGN sample is likely incom- 
plete because of spectroscopic incompleteness in the ChaMP 
and SEXSI surveys. The census of other cluster galaxies is 
very incomplete because few very high redshift clusters have 
the same quality membership data as our low-redshift sam- 
ple. In the first three subsections below we describe the choice 
of the fiducial absolute magnitude threshold, our estimate of 
the completeness of the spectroscopic observations of X-ray 
sources, and the total number of cluster galaxies in the clusters 
with incomplete membership data. The fourth subsection de- 
scribes our main result, the measurement of the AGN fraction 
and its evolution. The final two subsections describe potential 
contamination by AGN associated with large-scale structure 
around these clusters and other sources of uncertainty, respec- 
tively. 

4.1. Host galaxy magnitude threshold 

In previous work we defined the AGN fraction in clusters 
relative to galaxies more luminous than an j?-band ab solute 
magnitude of Mr = -20 mag (e.g. iMartini et ail 12006). This 
choice of magnitude threshold was largely driven by expedi- 
ence, namely it corresponded to the completeness limit for 
the most distant clusters in that sample. To properly extend 
this work to high redshift it is important to account for the 
evolution of the galaxy population in clusters, both in lumi- 
nosity and number. These were not significant effects in our 
low-redshift study as the highest-redshift clus ter was at only 
z = 0.3 1 , but in our previous work at z ~ 0.6 by lEastman et al.1 
d2007b the Mr = -20 mag cutoff corresponded to a fainter ab- 
solute magnitude relative to M^. Because the cluster galaxy 
population is larger, this would have led to a lower estimate of 
the AGN fraction if the cluster AGN are predominantly asso- 
ciated wit h the most luminou s galaxies, as is the case at low 
redshifts dSivakoff et al.ll2008l) . 

Here we adopt an absolute magnitude threshold of Mr(z) + 
1 , and thus allow for evolution of M%. At lo w-redshifts 
(0.01 < z < 0.07) IChristlein & Zabludoff (2003) measured 
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TABLE 6 
NONMEMBER X-RAY SOURCES 



CXOU ID 


z 


z ref 


R 


Rflag 


fx.s 


fx. Hi 


fx.B 


log L x ,b 


log L X .H 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


(10) 



Jl 12314.9+431208.3 


0.08017 ±0.00010 


1 


17.66 ±0.08 





Jl 12357.4+431314.1 


0.08007 


2 


19.46 ±0.08 





Jl 12403.0+431330.6 


1.1049 


2 


18.39±0.08 





Jl 12413. 1+430639.3 


2.3666 ±0.0015 


1 


19.80 ±0.08 





J 143804.9+033752.6 


0.29192 ±0.00030 


1 


18.50 ±0.06 





J143832.2+033506.0 


1.0083 ± 0.0051 


1 


19.98 ±0.06 





J143833.0+033606.8 


0.38252 ± 0.00017 


1 


19.40 ±0.06 





J143839.7+033631.3 


2.1493 ±0.0019 


1 


19.00 ±0.06 





J143841.9+034110.2 


1.7372 


2 


17.82 ±0.06 





J143847.3+032950.8 


0.00034 ±0.00012 


1 


16.89 ±0.06 





J 143859.0+033547.8 


0.7339 


2 


18.51 ±0.06 





J145623.0+221833.5 


0.00027 ±0.00010 


1 


15.51 ±0.07 





J145624.5+222057.1 


0.00019 ±0.00010 


1 


15.45 ±0.07 





J145634.6+221514.2 


0.40918 ±0.00010 


1 


20. 16 ±0.07 





J145657.7+221315.6 


0.00016 ±0.00010 


1 


14.87 ±0.07 





J145708.7+222352.4 


0.1238 


2 


17.44 ±0.07 





J145710.7+221844.9 


1.885 ±0.0014 


1 


18.73 ±0.07 





J145712.3+221446.7 


-0.00069 ±0.00010 


1 


15.15±0.07 


1 


J145721.0+222334.5 


1.7362 ±0.0010 


1 


19.33 ±0.07 


1 


J145726.9+221755.1 


1.4664 ±0.0011 


1 


19.55 ±0.07 





J151427.0+363803.1 


0.1616 


2 


16.90 ±0.06 





J151428.4+363743.5 


0.4026 


3 


20. 13 ±0.06 





J151437.5+364041.3 


0.1468 


3 


19.86 ±0.06 





J153938. 1+662102.4 


0.4375 


4 


19.71 ±0.08 





J154012.3+661439.2 


1.0577 ±0.0029 


1 


19.75 ±0.08 






8.4+1-5 
23.8+H 
22.2+I! 
7 53 +i:3 

J - z -0.95 

65 7+ 3 - 2 
OJ - -3.0 

11 5+ 1 - 3 

28.2+f 3 
16.3+*! 
46.5+1° 
9.<t« 
14.8+1* 
25.8+!^ 

■0.57 

50.4+2- 1 
9.4+' : A 



2.27 
3.99 



n 9O+0.6I 
2 - 2S -0.49 

12.1!!:, 

5 02+ 1 - 2 
J - u -0.98 

37.2+1? 



29 q+4.7 

17 3+*- 

7 56+ 2 5 
' - J -2.0 

<4.0 
79.8+f; 2 

1 8.3^'l 

16-l| 

29.5+| : 

2.0+ 2 - 2 

7 1+^ 

5 ^ 
3,1 -2.g 

56.8+« 
3.6+j;f 
< 3.4 

5 3 +L4 

15 3+ 2 - 2 
J ' -2.0 

7.0+" 
33.0+J ; | 

14. 

13.5_2 g 

6 5+ 2:< > 
41.9+" 



30. 1+ 3 ' 4 
55.9+^ 



++. 



44.5_ 4 
16.1+ 2 



+f 



149.3+f 
28.5+ 2 
34.8" 
61.6^ 
26.9+J 

1 13.8+1 
17.1+J 
29.8+? 
73.3 
16.0 



3 32+ 1 
9.4+[ 

90.2+ 3 

18.9+f 

56.3 
4.82 

18. 

26.9+| : f 

1 1 7+ 2:3 
—2.1 

83.1+1? 



■2.0 



41 69+ 005 

-Rftt 
44 51+0- 03 
-Q-P.3 

43 1 6+0"* 

44 73+0.07 

-8-Ofi 

4 2 32+0.10 

44 o<:+0.02 

~Q-Q? 

43 15+ 004 

-P.Q4 
44 .97+0.04 

45 01+^ : ^ 
44.41+0R2 



43 fi 3+0.02 
^ J - DJ> -0.02 



41.14+0 .13 
44 IB-t&Oo 
*+ ,zo -0.05 

44 <50+ 05 
-9Pi 

44 8 1 +0 03 
-Q.02 

41 53+0- fl 

4301*8 

^ JAn -0.05 

42 in+0.04 

42 on+0 09 

44.64+00? 



41 50+0 .05 
44 32+O.M 

42 go+0.04 

44 =54+0.07 
-Q06 
42 1 3+0-10 
44 66+0.02 
+^.°°-0.02 
42 Qfi+0-04 
-Q-Q4 
44 79+O.04 

4^83^1 



44.22+ 



43.45+ 



40 95+0 .13 

44 nQ+0- 06 
^• uy -0.05 

44 3 2 +o.05 
-Q.Q4 

44 62+0- 03 
-Q.02 

41 S5+0' T 

42.8311 
42.01+0 : 4 

1-0 09 
'-0.08 

44.46+ 002 



42.71+ 



"-0 07 



NO TE. — Chandra Observatio n L og. Columns are: Col (1) Name of X-ray so urce; Col (2) Redshift Col (3) References for redshift are: 1: this work; 2: 
SDSS I Adelman-McCarthy 2051) : 3: lAbraham et al] fT99"8l) ; 4: IMiller et alj 00041) ; Col (4) R-band magnitude; Col (5) Flags for photometry are: (0) no flag; 
(1) may be contaminated by nearby neighbors or bad pixels; Cols (6-8) Soft [0.5-2 keV], Hard [2-8 keV], and Broad [0.5-8 keV] band flux in the observed 



rest-frame in units of erg s 1 corrected for Galactic absorption. We do not quote luminosities for X-ray sources identified with Galactic stars (z ~ 0). 

the .R-band luminosity function (LF) for six nearby clusters 13 
and found that the composite cluster LF is consistent with a 



Schechter function with M* R = -2 1 .92 ± 0. 17 mag {h = 0.7, a = 
—1.21). They also find an essentially identical value of M R = 
-21.93 mag for the field. The low-redshift value of + 1 
is therefore about one magnitude brighter than the value of 
Mr = -20 mag we adop t ed in our previous, lo w-redshift stud- 



ies ( Martini et al.l 
son, Blanton et al. 



200a ISivakoff et alj|2008b . For compari- 
(120031) measured M* = -2 1 .22 (a = -1 .05) 
at z = 0. 1 for the r 01 band on the AB system. This corresponds 
to = -21.72 mag on the Vega system for the j?-band at 
Z= b ased on the conversions presented in Blanton & Roweis 
( 20071) and is therefore consistent with Christlein & Zabludoff 
d2003h . 

Many recent studies have measured the evolution of M R 
and generally these measurements include both a value for 
all galaxies and separate measurements for particular spectro- 
scopic types. This has relevance for our study as the clus- 
ter galaxy population is on average more quiescent than field 
galaxies and consequently their evolutionary history is differ- 
ent. We are most interested in measurements of the evolution 
of Mr as a function of spectral type to isolate the evolution 
of galaxies dominated by older stellar populations that are 
most likely representative of the evolution of cluster galaxies. 
A useful, low-redshift bench mark for a type-dependent LF 
for clusters comes again from IChristlein & Zabludoff (2003). 
They find M\ = -21.78 mag for quiescent galaxies in clus- 
ters, which is nearly identical to the value for all cluster mem- 

13 Tw o of these clusters ( Abell 85 and Abell 754) are in our low-redshift 
sample ( Sivakoff et alJ2008ft , 



bers. For field galaxies Che n et aTl d2003l) use photometric 
redshifts in the Las Campanas Infrared Survey and measure 
values of -21.70 to -22.22 mag (a = -1) for all galaxies over 
the range 0.5 < z < 1.5 and values of -21.21 to -21.82 mag 
(a = -0.2) for galaxies co n sisten t with an E/S0 + Sab spec- 
tral template. Wolf et ali (120031) use photometric redshifts 
from COMBO- 17 and measure more pronounced evolution 
for their early-type spectral template with M% fading by ~ 1 
mag from z~ 1.1 toz — 0.3. More recently. lllbert et~atj (2005) 
measure a fading of 1 . 1 - 1 .8 mag between z ~ 2 and z ~ 0. 1 in 
the 7?-band based on spectroscopic redshifts, although they do 
not present the evolution as a function of spectral type. These 
measurements of evolution in M\ are broadly comparable to 
the 1 .2 mag of fading from z = 1 to the present expected from 
pure luminosity evolution of a single stellar popul ation with 
Z/ = 2 and solar metallicity dBruzual & Charlol2003l) . 

Direct measurements of evolution of the clust er LF have 
mostly been conducted in the rest-frame B-band. iGoto et aTl 
(2005) find M* = -21.13 mag for MS1054-03 (z = 0.83), 
which is in our sample, and similar to the M% = -21.15 
mag measured for th ree clusters at an average z = 0.859 by 
Postman et al. (2001J. In compar ison to local fl-band mea- 
surements of the cluster LF (e.g. ICollessI fl989b iRauzy et all 
119981) . IGoto et al.l (120051) conclude that M* B fades by 0.46 to 
0.71 mag between z = 0.83 andz = 0. For t he same simple stel- 
lar po pulation model considered above dBruzual & Charloll 
2003), 1 .2 mag of fading in B-band is expected from z = 0.83 
to the present. While there is not a direct measurement 
in the rest-f rame i?-band for th e cluster LF, at yet longer 
wavelengths lEllis & Jonesl (120041) find that the fading in the 
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A'-band is 1.2 mag from z = 0.9 to the present and consistent 
with passive evolution and a formation epoch at z/ = 2. From 
these investigations of the LF evolution in the field and clus- 
ters, we adopt the assumpt ion that Mjjz) = M%(0)-z and the 
normalization for M% from Christlein & Zabludoff (20Q3|) for 
all cluster galaxies to estimate the completeness of the spec- 
troscopy of X-ray counterparts and the size of the galaxy pop- 
ulation in low-redshift clusters. This result is broadly consis- 
tent with all of the results described here, although is most 
consistent with the studies that predict more fading. If there 
is less fading of galaxies at the bright end of the LF, such as 
may be due to some low-level star formation in these galaxies, 
then the completeness limits we describe next are too bright 
and we will have systematically underestimated the popula- 
tion of luminous AGN in the higher-redshift clusters. 

4.2. Completeness 

We calculate a completeness limit in the observed /?-band 
for each cluster based on the value of M^(z) + 1 and a 
^-correction derived from the elliptical t emplate of the four - 
component spectral template presented bv lAssef et alj d2008l) . 
These templates are derived from 16,033 galaxies with spec- 
troscopic redshifts and multiband photometry from the AGN 
and Galaxy Evolution Survey. Most of the galaxies are in the 
range < z < 1 and the median redshift is 0.31. The par- 
ent sample is therefore broadly representative of our redshift 
range. For the higher-redshift clusters the ^-correction re- 
quires a substantial extrapolation from the observed /?-band, 
which for example samples rest-frame fi-band at z = 0.5. Our 
assumption that the typical cluster galaxies are best approxi- 
mated by an elliptical template is certainly reasonable for the 
low-redshift clusters. This may not be as good an approxi- 
mation at higher redshifts, although in a study of the color- 
magnitude relation in our tw o highest-redshift clusters (Lynx 
E and W) iMei et all d2009t) found there is no evidence for 
significant evolution. If a later-type template were a better 
choice for the ^-correction at higher redshift, the ^-correction 
would be smaller and the necessary /?-band spectroscopic 
limit would be brighter. The net effect would be a smaller 
completeness correction. 

The spectroscopic completeness of the high-z AGN sam- 
ple largely depends on the completeness of the ChaMP and 
SEXSI surveys, although we also use additional spectra for 
MS 2053.7-0449, MS 1054-03, and RDCS J0910+5422. The 
ChaMP survey quot es a spectroscopic c ompleteness of 77% 
for R< 22.37 mag dSilverman et al.l2005l) and the SEXSI sur- 
vey quotes a spectroscopic completeness of 61% for sources 
with 22 < /? < 23 mag, 67% for sources with 23 < R < 24, 
and 74% for source s with R > 24 mag (typically to 24.4 mag) 
(lEckartetal .1120061) . For the ChaMP data we adopt 77% as 
the completeness correction for 7? < 22.37 mag, while for the 
SEXSI survey we adopt an average completeness correction 
of 67% for R < 24.4 mag. For nearly all of the clusters above 
z > 0.6 the spectroscopic data do not extend to the equivalent 
of Mr(z)+ 1 and the size of the magnitude range without spec- 
tra ranges from a few tenths to over a magnitude. To estimate 
the number that may have been missed we inspected the host 
galaxy absolute magnitude distribution of the Lx.h > 10 43 erg 
s" 1 AGN in the clusters with complete data and find only one 
AGN fainter than M^. The distribution in M R of the X-ray 
AGN is shown in FigureQ] We therefore assume that we have 
not missed any AGN because the spectroscopic observations 
of X-ray sources did not have the requisite depth, although 
this assumption may have led us to underestimate the AGN 




-1 -2 -3 -4 
M R - (M* + 1) 

FIG. 1 . — Distribution in absolute magnitude Mr of the cluster AGN relative 
to M£ (z) + 1 at their redshift. All of the cluster AGN are substantially brighter 
than Mft (z) + 1 , although in most cases the spectroscopy is complete to this 
limit. The subset that are classified as BLAGN are represented by the hatched 
histogram. The dotted line corresponds to our galaxy luminosity threshold at 
M*(z)+\. 

fraction at high redshift. In contrast, if our assumption of an 
early-type template for the fe-correction was too red, then the 
spectroscopic data do achieve the requisite depth and this re- 
mains a nonissue. At brighter apparent magnitudes we do ap- 
ply a completeness correction to account for the quoted 77% 
and 67% completeness of the surveys. We discuss this further 
in §g3]below. 

The X-ray AGN populations of several of these clusters 
have been studied in previous work. The first substantial 
study of spectroscopic ally-confirmed X-ray AGN in a high- 
redshift cluster was by Uohnson et alj (T2003) in MS 1054-03. 
They identified 2 AGN associated with this cluster: CXOU 
J105702.7-033943 and CXOU J105710.6-033500; however, 
neither of these are included in the present sample. The 
first was not included because the X-ray luminosity is be- 
low our threshold of 10 43 erg s -1 and the second because it 
falls slightly outsid e the projected virial radius (R/R200 = 1 -2). 
iMartel et al.1 d2007l) have also studied X-ray sources in clus- 
ters, including three clusters that overlap this sample. They 
are discussed further in ^5.3| below. 

4.3. Inactive Cluster Galaxy Population 

To estimate the AGN fraction in these clusters we need 
to know the number of cluster galaxies more luminous than 
Mff(z)+1. We estimate this quantity in two ways, de- 
pending on the available data for the cluste rs. For the 
low-r e dshift clusters in our prev ious studies dMartini et al.l 
120061 120071: ISivakoff et all 120081) we have a large number 
of spectroscopically-confirmed cluster members and can es- 
timate the number of cluster galaxies either directly or with 
a completeness correction. We have calculated new estimates 
for these clusters for the present paper because we no longer 
use the Mr = -20 mag threshold of the previous work. These 
values are listed in Table [7] 

For essentially all of the new clusters in the present study 
we employ the same technique as East marTet al.l (2007) to es- 
timate the number of cluster members above M R > (z)+ 1 from 
the cluster velocity dispersion. This employs the richness- 
velocity dispersion relationship defined by iKoester et al.l 
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d2007l) for the MaxBCG cluster sample. The cluster rich- 
ness is the number of red (E/SO) cluster members more 
luminous than 0.4L* within the projected R200 radius. This 
relationship was originally derived from a sample of 13,823 
clusters with 0.1 < z < 0.3 in the SPS S with veloc i ty dis- 
persions greater than ~ 400 km s -1 . iBecker et al.l (2007) 
provide the most recent estimate of this relation based on 
a larger sample that extends over both a broader redshift 
range and to lower velocity dispersion groups. They find 
lncr = (6.17 ±0.04) + (0.436 ± 0.015) lnA^ a f°/25. For refer- 
ence a 520 km s" 1 cluster has N&f = 30. 

There are several caveats that need to be considered with 
the use of this estimator. First, the richness-velocity disper- 
sion relationship is based on photometric and not spectro- 
scopic redshifts. This is not a significant concern because 
for red cluster galaxies the photometric redshift estimates are 
robust within the quoted uncertainties. The second concern 
is that this relationship is based on the red cluster galaxies 
alone. At low redshifts this estimate is a reasonable approx- 
imation as the vast majority of cluster galaxies more lumi- 
nous than M* R + 1 fall in this category. For example, the frac- 
tion of quie scent galaxies above thi s lum inosity in the com- 
posite LF of IChristlein & Zabludof fl d2003l) is - 85%. While 
their definition of quiescence is based on spectral lines rather 
than color, these two definitions of quiescence typically agree 
when averaged over a cluster. At higher redshifts a larger 
fraction of the cluster galaxies may be blue due to ongo- 
ing star formation, but this can not be a substantial contri- 
bution because the luminosity-weighted mean star formation 
epoch is z = 2 for early-type cluste r galaxies up t o z = 0.5 
(Ivan Dokkum & van der Marel 2007). B ecker et al.ld2007l) do 
find evidence of evolution in this relationship in the sense 
of lower richness at fixed velocity dispersion in higher red- 
shift clusters, but they note that this may be due to their strict 
color selection. In addition, for our accounting of the inac- 
tive galaxy population the color of the galaxies does not mat- 
ter so long as they are in the cluster and above the luminos- 
ity threshold. Observations of individual clusters with exten- 
sive spectroscopic data support the assumption that there is 
no substantial evolution in the r elation between ha lo occupa- 
tion number and cluster mass (Muzzin et al. 2007). This is 
also supported by several theoretical studies that find minimal 
evolution in the num ber of bright galaxi es in massive halos 
dKravtsov et al.ll2004tlZentner et alj|2005l) . 

We performed an independent validation of the MaxBCG 
relation with an analysis of the individual clusters in our sam- 
ple with substantial membership data. While most of the low- 
redshift clusters have substantial membership dat a, these data 
gener ally do not extend to our estimate of R200 (Marti nTet al.l 
2007), nor is the X-ray coverage complete to this radius. Our 
spectroscopic coverage was often limited to the size of the 
Chandra field of view. However, two useful exceptions are 
Abell 89B and MS 1008. 1-1224 and in both cases estimates 
agree to within a factor of two. Our wide-field X-ray cov- 
erage of Abell 85 and Abell 754 (ISivakoffet al.|[2008h were 
designed to sample a substantial fraction of the projected R200 
and these values also agree well. Figure |2] illustrates the dif- 
ference between the MaxBCG membership estimates and our 
spectroscopic estimates. The larger points have nearly com- 
plete spectroscopic coverage to R200, while smaller points are 
substantially more incomplete. These points indicate that the 
error introduced by adopting the MaxBCG relation is approx- 
imately a factor of two. This error estimate is also consistent 
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FIG. 2. — Difference between predicted and measured cluster richness com- 
pared to the cluster richness predicted by the MaxBCG sample. The quantity 
is the number of red cluster galaxies more luminous than 0.4L* and es- 
timated from the cluster velocity dispersion 1 Beck er et al. 2007), while N spec 
is a spectroscopic estimate of this quantity (see j|4.3> . Symbols are coded 
according to the spectroscopic completeness relative to R200- Large circles 
have complete coverage to i?200> medium circles have more than 50% cov- 
erage, and the small circles have less than 50% coverage. Most clusters are 
at z < 0.5 (open sym bols), although substantial data exist for three at z > 0.5 
(filled symbols). See j|4.3| for further details. 

with an examination of figure 4 of Becker et ail (120071) . 

At higher redshifts three of our clusters have ex- 
tensive membership information. We estimate that 
MS 0015.9+1609 has ~ 200 members based on sev eral stud- 
ies dDressler & Gunnl Il992t lEllingson etaf] 119981) an d that 
MS2053.7-0449has~ 100 members dTran et al.l l2005). Note 
that these estimates are different from those presented in 
Eastm arTet al.l d2007l) due to updated completeness correc- 
tions and the change in the absolute magnitude threshold. 
For MS 1054-03 we estimate that there a re ~ 300 member s 
from the extensive spectroscopic work of iTran et al.l d2007l) . 
These three clusters are also shown in Figure[2](/;//<?ii circles). 
They are consistent with the low-redshift results and a factor 
of two uncertainty in the richness - velocity dispersion rela- 
tion. While our estimates of the cluster galaxy population for 
these three clusters, as for the low-redshift clusters, are based 
on all galaxies rather than just red galaxies, the consistency 
supports the assumption that the integral of the bright end of 
the galaxy luminosity function in clusters above an evolving 
Mr threshold scales reasonably well with the cluster velocity 
dispersion independent of redshift, even if there is evolution 
in the colors of the cluster galaxies. The number of AGN, 
estimate of the inactive population, AGN fraction, and spec- 
troscopic completeness for each cluster is listed in Table [7] 

4.4. Cluster AGN Fraction and Evolution 

The AGN fraction for any single cluster is very small and 
it is uncertain due to small number statistics. In addition, the 
AGN fraction may vary from cluster to cluster due to corre- 
lation s with other cluster properties such as velocity disper- 
sion dSivakoff et al .1120081) . The AGN fraction may also de- 
pend on variations in the properties of the galaxy population 
within each cluster (e.g., mass, SFR, morphology). We there- 
fore have binned the cluster sample in two ways to character- 
ize variations with redshift. First, we simply split the sample 
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FIG. 3. — Evolution of the AGN population in clusters from z = to 
z = 1.3 (filled symbols). The fraction of cluster members more luminous 
than M£ + 1 with AGN that have Lx.h > 10 43 erg s~'is shown in two red- 
shift bins (z < 0.4, z > 0.4; filled circles) and three redshift bins (z < 0.3, 
0.3 < z < 0.8,0.8 < z < 1.3; filled triangles). We also show our estimate 
of the field AGN fraction based on th e galaxy LF estimates by 
20051 open triangles), Dahlen et al. ( 20051 open circles), and 
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2003, open squares). See j|4.4| for further details. 



Chen et al 



at z = 0.4. This choice is primarily motivated by the tran- 
sition between where we rely on our own measurements and 
where we largely rely on other work. It also approximately di- 
vides the sample in two (17 clusters at z < 0.4, 15 at z > 0.4). 
This yields completeness-corrected AGN fractions of f A (z = 
0.19) = 0.00134!°°°^ and f A (z = 0.72) = 0.0100^;°°|, or ap- 
proximately a factor of eight increase in the AGN fraction (see 
Table [8) from a median redshift of 0.19 to a median redshift 
of 0.72. AGN fractions without the completeness correction 
are also listed in Table [8] The uncertainties on these quanti - 
ties are double-sided, 1 — a confidence limits (Gehrels 1986). 
The increase in the AGN fraction is formally significant at 
the 3.8cr level. We also split the sample into three bins with 
z < 0.3, 0.3 < z < 0.6, and z > 0.6 to better resolve the con- 
tinued increase at high redshift that is apparent in the raw data 
for individual clusters. This binning yields AGN fractions of 
f A (z = 0.15) = 0.0019!°;°°% f A (z = 0.45) = 0.003 1+™™, and 
f A (z = 0.81) = 0.0147!°;°°^. The measured evolution between 
the lowest and highest bins is also a factor of eight and in good 
agreement with the other binning scheme. We note that the 
observed evolution is also well fit by a simple power law scal- 
ing as f A oc (1 +z) a where a = 5.3+j although the power-law 
index is strongly correlated with the z = value of the AGN 
fraction. 

The factor of eight evolution of the AGN fraction is smaller 
but consi stent with the order of magnitude evolution ob- 
served by East man et al.l (120071) . They measured f A (z = 0.2) = 
0.0007:™ and f A (z = 0.6) = 0.020^ for L X , H > 10 43 erg 
s" 1 , although for a lower and fixed galaxy absolute magnitude 
of Mr = -20. At z = our galaxy absolute magnitude thresh- 
old is approximately a magnitude brighter than that used by 
Eastm an et al.1 (120071) and the offset increases linearly with 
redshift. This difference in absolute magnitude threshold can 
readily account for the change in the low-redshift fraction be- 
cause most of the AGN are associated with luminous clus- 
ter galaxies, that is increasing the galaxy luminosity threshold 
decreases the denominator and does not affect the numerator 
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2 
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1.00 
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15 


2 
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1.00 
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3C295 


0.460 


1642 


2 


412 


2 


n 4Q+0.64 


0.67 


MS045 1.6-0305 


0.538 


1371 





273 


2 


< 0.42 


0.77 


Moou i d . y + 1 buy 


0.541 


1234 


1 


214 


2 


n a -7+1.1 
0.47 


0.77 


KAJUo4o. /-HHf JO 


yj.j /*+ 




1 
J 


J uz 




A QO+2.3 

0.98 8[ 


u.o / 


MS2053.7-0449 


0.583 


865 





95 


2 


< 1.2 


1.00 


RXJ0542.8-4100 


0.634 


1269 


5 


229 


2 


2 - 18 -oii 


0.77 


RXJ2302.8+0844 


0.722 


658 





50 


2 


< 2.3 


0.77 


MSI 137.5+6625 


0.782 


885 


1 


100 


2 


l.oo^l 


0.77 


RX J1317.4+2911 


0.805 


1142 


1 


179 


2 




0.67 


RXJ1716.4+6708 


0.813 
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1 
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NOTE. — AGN fraction estimates for individual clusters. Columns are: Col. (1): Cluster 
name; Col. (2): Redshift; Col . (3): Velocity dispersion (references for these values are in Tablefn 
Table[3l [Sivakoff et alj(200l for Abell 754, Abell 85, Abell 89B, lMartini etlul t2Tjfj7J) for Abell 
3128, Abell 3 125, Abell 644, Abe ll 2104, Abell 2163 and M S1008.1-1224, or adopted from 
Czoske i 2004J) for Abell 1689 and lGirardi & Mezzett) t200il) for AC 114); Col. (4): Number 
of AGN with Lx h > 10 43 erg s~4n galaxies more luminous than M^(z)+ 1; Col. (5): Estimate 
of the number of cluster galaxies more luminous than (z) + 1 within either the Chandra FOV 
or R200, whichever is smaller; Col. (6): Flag for the origin of the estimate where 1: from our 
spectroscopy and completeness correction; 2: from the MaxBCG as described in j j4.3| Col. 
(7): Estimate of the cluster AGN fraction in percent; Col. (8): Estimate of the spectroscopic 
completeness for X-ray sources. 

of the AGN fraction. In addition, we h ave since identified a 
second luminous AGN at low redshift (Sivakoff et al. 2008). 
At high redshift the change in galaxy luminosity threshold is 
also important, but in addi tion the cluster sample is more than 
three times larger than the lEastman et alj d2007l) sample. The 
low-redshift cluster sample has increased by less than a factor 
of two. 

One way to characterize the evolution of the cluster AGN 
fraction relative to the field is to calculate the integral of 
the field space density $(L^.// > 10 43 ) as a function of red- 
shift. Integr ation of the l umino sity-dependent density evolu- 
tion model in lUeda et alj {2003) yields a factor of five increase 
between z = 0.8 and z = 0.2, which is somewhat less but con- 
sistent with the observed evolution of cluster AGN. However, 
this is not a fair comparison because the evolution of field 
AGN with ${Lx.h > 10 43 ) is not normalized by the evolution 
of all field galaxies brighter than Mr + 1 and the cluster AGN 
fraction is. 

While there is not a direct measurement of the field AGN 
fraction similar to our calculation for clusters (although see 
Lehm eret alj2007h . we can estimate this qu antity by d i viding 
the integral of the field hard X-ray LF from lUeda etal] d2003l) 
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TABLE 8 

AGN Fraction for Subsamples of the Clusters 



Sample 
(1) 


Z range 
(2) 


Net 
(3) 


median z 
(4) 


median a 
(5) 


N A ,rm 

(6) 


Ngal 
(7) 


fA,raw [%] 

(8) 


fspec 

(9) 


/a.coit [%] 


Two Bins 


z<0.4 
Z > 0.4 


17 
15 


0.19 
0.72 


993 
895 


2 
18 


1492 
2379 


134+0- 18 
u - 10 -0.087 

76 +0 22 

u ' -0.18 


1.00 
0.76 


M1.087 
1 00+ ' 29 


Three Bins 


z<0.3 

0.3 <z<0.6 

z>0.6 


13 
10 

9 


0.15 
0.45 
0.81 


953 
1065 
885 


2 
4 
14 


1049 
1604 
1218 


qj 0+0.25 

25+^ 


1.00 
0.81 
0.78 


19+0.25 

0.31^41 

1 47+8:50 



NOTE. — Cluster AGN fractions with the data split into two bins and three bins. The two bins are split at z = 0.4, while the three bins split the data at z = 0.3 and z = 0.6. For each 
bin we list: Col. (2): redshift range; Col. (3): number of clusters; Col. (4): median redshift; Col. (5): median velocity dispersion of clusters; Col. (6): sum of the luminous AGN in 
the bin; Col. (7): raw AGN fraction with double-sided, 1 - a confidence limits; Col. (8): estimate of the mean spectroscopic completeness weighted by the number of galaxies per 
cluster; Col. (9): AGN fraction corrected for spectroscopic completeness. 



by the integral of the galaxy LF. We have identified three sur- 
veys that report LF measurements for the /?-band and approx- 
imately span the same redshift range of t his work. The firs t 
of these is the VIMOS-VLT Deep Survey (lllbert et alJl2005h . 
which is based on UBVRI photometry, ~ 1 1 , 000 spectra to 
Iab = 24 mag and extends from z = 0.05 to z = 2 (although 
their lowest-redshift point is taken from SDSS: lBlanton et all 
2003). We also show results from two measurements based 
on photometr ic redshift data: th e Las Campanas Infrared Sur- 
vey (LCIRS: IChen et alj I20031) . which is mostly based on 
U BVRIH measurements and presents the LF for z = 0.5— 1 .5, 
and the Great Obse rvatories Origins Deep Survey (GOODS; 
iDahlen et alj|2005l) . which is based on U through K obser- 
vations and presents the galaxy LF to z = 2. While these 
photometric redshift surveys may have more systematic un- 
certainties than the LF based on spectroscopic measurements, 
they have the virtue that they have measured the luminosity 
function in the rest-frame /?-band, rather than relied on as- 
sumptions about galaxy spectral energy distributions (SEDs) 
to calculate ^-corrections. We have calculated the field AGN 
fraction for each of these surveys and show the results in 
Figure [3] {open sym bols). At low reds hift the AGN fraction 
calculated with the lllbert et all d2005l) LF is approximately 
a factor of five above the cluster fraction, which is consis- 
tent with the d i fferen ce between the field and clusters seen by 
iDressler et alj (1999) for spectroscopically-identified AGN. 
At higher redshifts (z > 0.5), the field estimates range be- 
tween a factor of three and a factor of ten above the cluster 
fraction. These estimates of the field AGN fraction vary so 
substantially due to the dispersion in estimates of the galaxy 
luminosity function. In addition, this calculation presupposes 
that all of the X-ray AGN are in galaxies more luminous than 
Mft(z)+ 1. While there is good evidence that most of these 
luminous X-ray AGN a re in relatively luminous galaxies (e.g. 
ISilverman et alj|2009al) . there is nevertheless a bias against 
spectroscopic identification of lower-luminosity X-ray AGN 
host galaxies. Finally, we note that the relative evolution of 
galaxies in clusters and the field further complicates this com- 
parison. In future work we hope to compile sufficient data 
to calculate the AGN fraction in the field and clusters as a 
function of galaxy mass. At present the data are insufficient 
to conclude if the cluster AGN fraction or field AGN fraction 
evolves more rapidly. 



4.5. Contamination by AGN Associated with Large-Scale 
Structure 



One concern raised about the physical origin of the 
Butch er-Oemler effect is th e contribution of projection ef- 
fects. iDiaferio et al.l (1200 ll) studied this issue in detail with 
N-body simulations and semianalytic models to distinguish 
true cluster members from field interlopers that were at the 
cluster redshift and with in the projected j?200 , yet physically 
outside the cluster /?200- Diaf erio et alj (J2001) concluded that 
up to 50% of the apparent Butcher-Oemler galaxies at the red- 
shifts of high-redshift clusters may be interlopers. A similar 
effect may be relevant for the AGN population and such a 
large contamination would decrease the observed evolution, 
but not erase it. 

While there is no comparable study that directly investi- 
gates the projection of AGN onto high-redshift clusters, there 
is good evidence that AG N are associated with t he large-scale 
environment of clusters. iGilmour et al.l d2007l) identified 1 1 
X-ray AGN (to a lower luminosity limit of ~ 10 41 erg s" 1 ) in 
the A901/2 supercluster at z ~ 0.17 and only one was in the 
densest region of the supercluster. The remainder were mainly 
in region s of intermed i ate de nsity. In the vicinity of 3C295 
(z = 0.46^ lD'Eliaet all (l2008l) find evidence for AGN associ- 
ated with a filamentar y structure. At yet hig her redshifts this 
trend is also apparent. iKocevski et a l. (2009) find X-ray AGN 
associated with the CL1604 supercluster at z ~ 0.9, which 
contains 8 confirmed groups and clusters. These AGN mostly 
avoid the densest regions of the clusters and are located on the 
outskirts of the most massive clusters, that is they are associ- 
ated with poorer clusters and groups. 

We examined our data to determine if there were a pop- 
ulation of AGN outside the projected 7?2oo for these clusters 
similar to those seen in the two superclusters. This is only 
possible with the subset of the sample with substantial cover- 
age beyond Z?200- Eight of the clusters have Chandra cover- 
age that extends to 27?2oo- There are six AGN between /?20o 
and 2/?2oo that meet our velocity cuts for cluster membership 
compared to eight AGN within /?200 for these same clusters. 
The larger number within the clusters suggests the opposite 
trend from the two supercluster studies described above, al- 
though these results are not truly in conflict because the su- 
percluster studies encompassed a much larger area outside of 
dense clusters than this study. The different large-scale envi- 
ronments associated with these clusters and the superclusters 
suggest a more quantitative comparison would not be mean- 
ingful. These large-scale structure data also provide a crude 
means to estimate the likelihood of chance juxtapositions of 
AGN associated with large-scale structure onto the clusters. 
If interloper AGN have the same surface density within /?200 
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as between /?2oo and 2/?200, then the six we identified in an 
area of 37i7?200 suggest we should expect at most 2 interlopers 
compared to the 8 AGN we see within /?200- This line of ar- 
gument suggests that the interloper fraction is 25%, which is 
small compared to the observed evolution signature. 

4.6. Uncertainties 

One major potential source of systematic error is the use 
of the MaxBCG richness estimator to estimate the fraction of 
cluster galaxies more luminous than + 1. In jj4.3| we esti- 
mated that there is a factor of two uncertainty in the use of this 
relation. This uncertainty is mainly important for the high- 
redshift subsamples as the low-redshift subsamples have more 
complete spectroscopic membership data. If we randomly in- 
troduce a factor of two uncertainty in each cluster, the effect 
is negligible when averaged over the 15 clusters with z > 0.4 
compared to the factor of eight evolution in the AGN fraction. 

As mentioned previously, another valid concern with the 
MaxBCG estimator is that it is calibrated to the number of 
red galaxies in the cluster and this population may not all be 
in place at z = 0.4 and higher. For our application it does not 
matter if the galaxies are red or not, just that they are in the 
cluster. Furthermore, if we have overestimated the number 
of galaxies brighter than M| + 1 then we have underestimated 
the evolution of the AGN fraction and our result is yet more 
statistically significant. The assumption that all of the galax- 
ies are red does impact the ^-correction we use to estimate 
the spectroscopic limit corresponding to M^(z)+ 1 and thus 
the size of our completeness correction. If the galaxies are 
redder, then the ^-correction would be smaller, the apparent 
magnitude limit would be brighter, and the completeness cor- 
rection would be smaller. The implication would be that we 
have preferentially overestimated the AGN fraction at high 
redshifts because completeness corrections are only applied 
to the high-redshift clusters. While the average completeness 
correction approaches 25% (see Table|8]l, in practice the spec- 
troscopic completene ss is not a strong function of apparent 
magn itude (e.g. see 34.21 Silverman et al. 2005; Eck art et"ai1 
2006) and we consequently expect much less than a 25% re- 
duction in the evolution. The evolution of the host galaxy 
population is also important because if there were less fading 
of M*(z) than we assume, then the completeness limit would 
be too bright and we would have underestimated the AGN 
fraction at high redshift. 

The value of the cluster velocity dispersion introduces addi- 
tional uncertainty to this calculation in two ways. First, many 
of the direct measurements of the cluster velocity dispersion, 
particularly for high-redshift clusters, are based on small sam- 
ples of galaxies and thus the velocity dispersion itself may 
be uncertain, particularly if the galaxy velocity distribution 
is not Gaussian. Second, as noted above the cluster velocity 
dispersion has not been directly measured for several clus- 
te rs and we instead u sed the X-ray temperature and the results 
of IXue & Wul (|2000) to estimate the velocity dispersion and 
this has a 30% scatter. We checked both of these concerns 
with a measurement of the scatter between a and Tx for the 
ten high-redshift clusters with measurements of both quanti- 
ties and the mean deviation is ~ 220 km s" 1 if we exclude 
3C295, whi ch has a substantially h igher velocity dispersion 
(1642 km/s lGirardi & Mezzettill200ll) than expected f rom its 
X-ray temperature (5.3 K from Vikhli nin et al.l 120021) . This 
mean deviation corresponds to approximately a factor of two 
uncertainty in the richness, which is comparable to the uncer- 
tainty we derived for the richness estimator. From this analy- 
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FIG. 4. — Histograms of the number of clusters with a given velocity 
dispersion (dotted line) and the number of AGN in clusters of a given veloc- 
ity dispersion (dashed line) for the low-redshift (z < 0.4; top panel) and the 
high-redshift (z > 0.4; bottom panel) subsamples. The cluster samples are 
reasonably well matched within these two redshift bins. 

sis we similarly conclude that this source of uncertainty does 
not substantially affect our results. 

A related evolutionary effect is that the velocity distribu- 
tions of the high-redshift clusters may be systematically more 
non-Gaussian than low-redshift clusters because the high- 
redshift clusters are less likely to be relaxed. If the cluster 
velocity dispersion were overestimated, then the richness and 
R200 would be overestimated as well. This in turn would lead 
to a n underestimat e of the AGN fraction in high-redshift clus- 
ters. iJeltema et al.l (12005) measured power ratios from Chan- 
dra observations of the IGM for a large sample of clusters 
out to z ~ 1 and found good evidence that high-redshift clus- 
ters are less relaxed than low-redshift clusters, so this poten- 
tial source of systematic error would lead us to underestimate 
the AGN fraction. Nin e of our clusters were analyzed in 
the IJeltema et al.1 (120051) study, including eight in our high- 
redshift sample. We compared the AGN fractions and the 
power ratios for these clusters, but did not find a significant 
trend. Unfortunately we do not have sufficient redshift data 
for most high-redshift clusters to look for non-Gaussianity in 
the galaxy velocity distribution, although note there is no evi- 
dence for a trend between dyn amically-disturbed clusters and 
AGN fraction at low redshift dMartini et al.ll2007t) . 

Finally, we consider the evolution of the cluster population 
to determine if the higher-redshift clusters represent the pro- 
genitor population of the lower redshift clusters. As noted 
previously, observations at low redshift indicate that the AGN 
fraction depends on environment and specifically that the 
AGN f raction is higher in lower velocity disp ersion environ- 
ments ( Siva koff et al1l2008l lArnold et al.ll2009h . Therefore if 
our high-redshift clusters are the progenitors of lower velocity 
dispersion clusters or massive groups, then the observed evo- 
lution may not be as significant. As many of the high-redshift 
clusters are X-ray selected, they are generally high-mass clus- 
ters and are reasonably well matched to the low er-redshift 
sample (se e Figure [H and Table |7J. Following iFinn et all 
(120051) and Poggianti et al.1 (120061) . we have estimated the ve- 
locity dispersions of the progenitors of the high-redshift clus- 
ter sample and find they are in good agreement. For exam- 
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TABLE 9 

High-Redshift AGN Associated with Large-Scale Structure around Clusters 



AGN 


Cluster 


z 


R [mag] 


logLx,// [ergs '] 


Sv/cr [km/s] 


Ai? [arcmin] 




Class 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(V) 


(8) 


(9) 


CXOSEXSI J084846.0+445945 


RX J0848.7+4456 


0.567 


21.45 


43.1 


1.99 


3.51 


1.16 


ELG 


CXOMP J230300.9+084659 


RXJ2302.8+0844 


0738 


21.71 


44.23 


2.81 


4.46 


1.2 


BLAGN 


CXOSEXSI J171 807.6+670647 


RXJ1716.4+6708 


0797 


21.75 


44 


1.83 


7.8 


1.59 


BLAGN 


CXOU J105710.6-033500 


MS 1054-03 


0.832 


21.93 


43.14 


1.27 


4.57 


1.18 


ALG 


CXOSEXSI J091040.8+542006 


RDCS J0910+5422 


1.097 


22.38 


43.1 


2.74 


2 


1.13 


ELG 


CXOSEXSI J084903.9+445023 


LynxE 


1.276 


23.92 


43.2 


2.95 


1.76 


1.03 


ELG 



NOTE. — AGN associated with large-scale structure around the subset of high-redshift clusters with complete X-ray coverage to twice the projected virial radius. This is the 
subset of AGN that satisfy the redshift selection criterion, but have a projected distance of 1 < R/R20Q < 2. Columns are id entical to Tabled The data for CXOU J 1057 10.6-033500 
are from van Dokkum et al. 1 2000) for the redshift, magnitude, and classification and the X-ray data are from Johnson et al. i 2003). This sample is described in further detail in 34.51 



pie, the progenitor o f a 1000 km s clu ster at z = has 800 
km s" 1 at z = 0.6 dPoggianti et al.ll2006l) . or only about 100 
km s -1 less than the difference between our low-redshift and 
high-redshift subsamples. The sense of this trend is that the 
high-redshift sample is actually somewhat more massive than 
the typical progenitor of the low-redshift sample and therefore 
the minor mismatch in cluster masses is more likely to have 
dampened rather than enhanced the measured evolution of the 
AGN fraction. 

5. PROPERTIES OF THE CLUSTER AGN 
5.1. Distribution 

The projected radial and velocity distributions of the AGN 
provide valuable additional information about the origin of 
the AGN. For example, if the AGN are preferentially located 
in the cluster outskirts, or preferentially have a higher velocity 
dispersion than the cluster mean, this may indicate that their 
host galaxies have relatively recently entered the cluster po- 
tential. This is known to be the case for e mission-line galaxies 
dBiviano et all 119971: iDressler et aD Il999h. At low-redshifts 
and for lower-luminosity X-ray AGN. iMartini et al.1 d2007l) 
found that L x > 10 42 erg s" 1 [0.5-8 keV] AGN were more 
centrally concentrated than typical cluster galaxies, while 
AGN an order of magnitude less luminous had the same distri- 
bution as the inactive galaxy population. For both luminosity 
thresholds the velocity distribution of the AGN were consis- 
tent with the galaxy population. 

It is more challenging to compare these higher-luminosity 
X-ray AGN to the host galaxy population because we lack 
membership data for nearly all of the high-redshift clusters. 
Nevertheless, we can compare the distribution of sources to 
the typical distribution of cluster galaxies and to the excess 
surface density distribution found by surveys of X-ray point 
sources toward distant clusters. In Figure we present a his- 
togram of the number of X-ray AGN from the cluster center as 
a function of distance in both physical units (Mpc) and nor- 
malized to 7?2oo- While the sample is small, two results are 
apparent from the figure. First, there are approximately equal 
numbers of AGN outside 0.5/? 200 as inside it, whereas if the 
AGN traced the cluster galaxy distribution we would expect 
them to be more centrally concentrated. Second, the radial 
distribution is more strongly peaked when plotted in physical 
units than normalized to R200; 

While we do not have detailed information on the radial 
distribution of the cluster galaxy populations in these clus- 
ters, we do have extensive data on nearby clusters from 
IChristlein & Zab ludoff (2003). For these clusters we have in- 
vestigated the cluster galaxy distribution with the same se- 
lection criteria (R < R 2 oo, M R < M| + 1, Av < 3a) and find 
that 70% of the galaxies fall within 0.5/?200, whereas 10 of 



18 luminous AGN at z > 0.4 are within 0.5/?200- The bino- 
mial probability is only 14% that we would find 10 or fewer 
AGN within 0.5/?200 if we expected 70%. There is thus a 
mild tendency for luminous AGN to be distributed toward the 
outskirts of the clusters, although this does make the substan- 
tial assumption that the radial distribution of galaxies within 
clusters is similar at z ~ 0.8 and the present. This broad 
distribution in radius is in contrast to our earlier results on 
lower-luminosity AGN in lower-redshift clusters. At low red- 
shift we found that 50% of the luminous AGN were within 
0.1/? 200 (Martini et a l.l2007l) . Better statistics could determine 
if the AGN are preferentially located in the outskirts of clus- 
ters compared to all cluster galaxies. That would be consistent 
with the hypothesis that A GN are trigg e red by mergers during 
infall. From simulations Ghig na et al.l ([1998) find that merg- 
ers betwe en galaxies d o not o ccur within the virial radius. We 
note that Berri eret al.1 (120091) simulated the formation of 53 
galaxy clusters and find most cluster galaxies do not expe- 
rience 'preprocessing' in group environments and therefore 
processes specific to clusters must largely be responsible for 
the differences between cluster and field galaxies. 

The second result has interesting implications for stud- 
ies that use the surface density distribution of ex- 
cess sourc es to characterize t he di st ribution of AGN in 
clusters dRuder man & E befingj 120051; iGilmour et all 120091; 
iGalametz et al.l 120091) . These studies generally plot the ex- 
cess surface density as a function of physical distance from 
the cluster center and find a central peak in surface density. 
Our results indicate that the true distribution may be flatter 
than implied by use of the physical (proper) distance from the 
cluster core. This is because those surveys, like the present 
study, include clusters with a wide range of masses and conse- 
quently a wide range of /?2oo- Simply adding the distributions 
for all clusters without renormalizing each observation for the 
size of the cluster will produce an artificial central peak due 
to the mass range of the cluster sample. 

If the cluster AGN are associated with a population that re- 
cently entered the cluster potential, the host galaxies may also 
be preferentially on more radial orbits and have a larger ve- 
locity dispersion than that of all cluster galaxies. As noted 
previously, this is true of the emission-line galaxy population 
in clusters. In Figure [6] we plot the cumulative velocity dis- 
tribution for all 18 AGN with z > 0.4 normalized by the clus- 
ter velocity dispersion. The distribution is in excellent agree- 
ment with a Gaussian distribution and we therefore find no 
evidence that the cluster AGN have a larger velocity distri- 
bution that would be consistent with more radial orbits. This 
was also found for the 14 relatively luminous (Lx.b > 10 42 erg 
s" 1 ) AGN studied bv lMartini et al.1 d2007l) . A better test would 
be to compare the AGN host population to the absorption-line 
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FIG. 5. — Histograms of the AGN clustercentric distances in terms of Mpc (left) and normalized to S200 (right) for cluster AGN with z > 0.4. The distribution 
of the confirmed cluster members (solid line) is much more centrally peaked when expressed in terms of Mpc than in terms of r/R^oo- Other AGN associated 
with large-scale structure (with R > R200) are also shown (dotted line). 
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FIG. 6. — Histogram of the cumulative velocity distribution of cluster 
AGN normalized to the cluster velocity dispersion for the 18 cluster AGN 
with 1 > 0.4 (solid histogram). The AGN velocity distribution is consistent 
with a Gaussian distribut ion (solid curve ) and the L x .b > 10 42 erg s _1 AGN 
fromfMartini et al. ( 2007) (dotted histogram). 



galaxies in the clusters since the velocity dispersion estimates 
for many of these clusters may be biased toward the emission- 
line galaxy population because it is easier to measure redshifts 
for them. While this is not the case for those whose veloc- 
ity dispersions are estimated from X-ray data, it may also be 
true of the calibration sample for the relations between X-ray 
properties and galaxy velocity dispersion. 

5.2. Luminosity Function 

We have begun to acquire sufficient numbers of cluster 
AGN that it is possible to compare the X-ray luminosity func- 
tion (XLF) between clusters and the field, as well as the clus- 
ter XLF at different redshifts. A comparison between the 
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FIG. 7. — Hard X-ray lu minosity function o f cluster AGN at z > 0.4 com- 
pared to the field XLF from Ueda et al. 1 2003) at the median cluster redshift 
(z = 0.8, solid curve) and at low redshift (z. = 0.1, dotted curve). The field 
XLFs have been renormalized to be consistent with the cluster measurements 
in the first two luminosity bins. The arrows are upper limits calculated with 
Poisson statistics. 

cluster and field XLF is interesting because differences be- 
tween the two would be a signature of environment-dependent 
downsizing. There is evidence that this is tr ue of star forma- 
tion in different environments. For example Kauffman n et alJ 
(2004) find that substantial star formation is only present in 
higher mass galaxies in lower density environments in the 
local universe. If the cluster black holes primarily grew at 
higher redshifts than field black holes, similar to the earlier 
formation epoch expected for the stellar populations in lumi- 
nous cluster galaxies, then the cluster luminosity function at 
high-redshift may have a similar shape to the present-day XLF 
in the field. One test of this hypothesis is to compare the char- 
acteristic luminosity L\ between clusters and the field. If the 
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cluster AGN primarily grew at an earlier epoch, L£ would be 
smaller in clusters relative to the field at a given redshift. 

It is reasonable to compare the shape of the XLF between 
clusters and the field because the XLF is a measurement of the 
X-ray sources alone within well-defined volumes, although 
t he c aveats associated with large scale structure discuss ed in 
jj4.5l do apply. This is different from the case in Section l4~4l 
where we noted that the comparison of the evolution of the 
AGN fraction and the integrated space density was not com- 
paring identical quantities because the AGN fraction includes 
information about the galaxy population. The one assumption 
that we do make is that all of the X-ray sources are hosted 
by galaxies above our threshold, but this is reasonable given 
FigureQ] In addition, the normalization remains arbitrary be- 
cause it is challenging to define a total volume for the cluster 
AGN sample, although this is not necessary because the shape 
of the XLF already provides useful information. In Figure [7] 
we plot the cluster XLF for our z > 0.4 sample compared to 
the field XLF at the median cluster redshift of z = 0.8 from 
lUeda et all (120031) . The cluster XLF is in reasonable agree- 
ment with the field XLF at the same redshift, although the 
statistics are quite limited. As motivation for future work, we 
also plot the field XLF at lower redshift (z = 0. 1, dotted line). 
For the lower-redshift XLF L\ is smaller and consequently 
all Lxfi > 10 43 AGN are above the characteristic luminos- 
ity, while these data straddle L\ in the field XLF at z = 0.8. 
Improved statistics for cluster X-ray AGN at z > 0.4 could 
determine if there is also a break in the cluster XLF, or if it is 
more similar to the field XLF at lower redshift. 

The evolution of the cluster XLF with redshift is also rele- 
vant for the origin of X-ray AGN in lower-redshift clusters. If 
cluster AGN at the present day are simply the descendants of 
AGN at higher redshift that have been fading for several Gyr, 
then the difference between the low -redshift and high-redshift 
cluster XLF should be consistent with pure luminosity evolu- 
tion. In contrast, if there is substantial retriggering of low- 
luminosity AGN in low-redshift clusters, or if other mecha- 
nisms are capable of fueling AGN in clusters, then the clus- 
ter XLF evolution may not be consistent with just luminosity 
evolution. A signature of other fueling or triggering mech- 
anisms would be a substantially larger population of lower- 
luminosity AGN in present day clusters compared to expecta- 
tions from the high-redshift population. While pure luminos- 
ity evolution would be surprising because this is not observed 
in field AGN, the most luminous cluster galaxies are consis- 
tent with passive evolution. Better measurements of the clus- 
ter XLF over a broader range in luminosity could investigate 
this hypothesis. 

5.3. Host Galaxy Properties 

Both the colors and morphologies of low-luminosity (~ 
10 41 erg s" 1 ) AGN in low-redshift clusters suggest they 
are primarily hosted by gala xies dominated by light from 
their old stellar populations (Mar tini et al.l 12006). This be- 
comes progressively less true for higher-luminosity AGN 
and ground-based observation s of the most lumino us sources 
(> 1 43 erg/s) in Abell 2 104 dMartini et al.| [2002) and Abell 
754 (lArnold et all |2009) indicate that they have late-type 
morphologies, although their hosts are luminous. In addi- 
tion, these more luminous AGN are more likely to exhibit 
visible-wavelength AGN spectral signatures than their lower- 
luminosity counterparts. 

While the spectroscopic classification of the high-redshift 
sample is fairly subjective because of variations in wavelength 



coverage and signal to noise ratio, the spec trosc opic classifi- 
cation s reported by ISilverman et al.l (120051) and lEckart et aD 
(2006) support the low-redshift results. They classified six of 
17 X-ray AGN as BLAGN, nine as other emission-line galax- 
ies, and the remaining two as absorption-line galaxies. The 
vast majority of the higher luminosity AGN have substantial 
line emission, even with the bias against redshift measure- 
ments for sources without strong emission lines. We note for 
comparison that two of the six AGN in the large-scale struc- 
ture sample are classified as BLAGN and the other four are 
evenly split between emission-line and absorption-line galax- 
ies. These other sources are thus similar to the cluster AGN. 

Several of the high-redshift clusters also have HST ob- 
servations suitable to study the morphologies of the cluster 
galaxies. The largest survey of X-ray source mor phology in 
high-redshift clusters is that by iMartel et al.l (120071) . who in- 
vestigate the fields of five high-redshift clusters: RX J0152- 
1357, RX J0849+4452, RDCS J09 10+5422, MS 1054-0321, 
and RDCS J 1252-2927, and the middle three clusters over- 
lap this sample. For the entire field sample they classify half 
of the X-ray counterparts as early-type, 35% as late-type, 
and 15% as irregular galaxies. For the six cluster members 
in their sample, they find half are in early-type hosts, two 
in late-type hosts, and one in an irregular galaxy. In ad- 
dition, three of these cluster AGN hosts are in interacting 
systems. The specific overlap with this sample are CXOU 
J09 1043. 3+542 152 and CXOMP J105650.6-033508 and both 
have early-type morphologies (their other member in RDCS 
J0910+5422 falls slightly below our luminosity threshold). 

5.4. Implications for the Sunyaev-Zel'dovich Effect 

Many cluster surveys are currently planned or in progress 
that use the Sunyaev -Zel'dovich effect to identify large num- 
bers of clusters (e.g. lKosowskyll2003l: iRuhl et alll2004l) . This 
effect is caused by inverse Compton scattering of Cosmic Mi- 
crowave Background (CMB) photons off hot electrons in the 
ICM that change s the spectrum of the CMB in the direction 
of a cluster (e.g. ICarlstrom et al.1 1200 2). The main virtue of 
this effect is that it is redshift independent, and consequently 
can be used to detect (the hot electrons associated with) clus- 
ters out to high redshifts. However, mechanical heating by 
AGN in the cluster may contr ibute to the thermal energy of the 
ICM (e.g. Birzan et al. 2004) and thus make it more difficult 
to identify some clusters. Any increase in the AGN popula- 
tion with redshift will also introduce a systematic effect with 
redshift. 

The potential impact of AGN on SZE cluster sur veys was 
recently examined in detail by iLin & Mohrl (120071) . They 
measured the radio luminosity function in nearby clusters at 
1 .4 GHz and used measurements of AGN at higher frequen- 
cies dCoorav et al.ll998l:lc"oble et al.l2007l) to estimate that on 
order 10% of clusters will have an AGN flux comparable to 
the SZE flux. As a worst-case scenario they adopted an evo- 
lution model where the fraction of radio AGN increases as 
(1 +z) 2 ' 5 - This model was largely motivated by observations 
of the radio galaxy lu minosity function, which suggested evi - 
dence for an increase (Be st et al.l2002tlBranchesi e t al. 2006). 
If this population evolves at a comparable rate more consistent 
with the (1 +z) 5 3 rate we observe for luminous X-ray AGN, 
then the fraction of substa ntially contaminate d clusters will be 
higher than predicted bv lLin & Mohrl d2007l) . 

6. DISCUSSION 
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The extent of the correlation between the evolution of star 
formation and AGN in clusters could provide valuable new 
insights into ho w closely related these two proces ses are. The 
original work bv lButcher & Oemlerfdl978l fl984) on the evo- 
lution of the fraction of blue galaxies in clusters provides a 
useful first point of comparison to the AGN fraction evolution, 
in part because we adopted many elem e nts of their method- 
ology. Specifically, Butc her & Oemleri (1 19841) characterized 
cluster galaxy evolution with: 1) a fixed criterion to define 
the sample of interest (a galaxy was classified as blue if the 
rest-frame B — V color was at least 0.2 mag bluer than the 
relation exhibited by the red galaxies); 2) measurement of 
this population relative only to cluster galaxies above some 
luminosity threshold (My = -20); 3) use of an aperture scaled 
to the physical properties of individual clusters (a circle that 
contained the inner 30% of the cluster g a laxy p opulation). 
With these definitions, iButcher & Oemleri (119841) found that 
the blue galaxy fraction increased from fs ~ 0.03 at z < 0.1 
to fs ~ 0.25 at z = 0.5 for relatively compact, concentrated 
clusters, or approximately an order of magnitude. 

One of the most recent and comprehensive studies of 
the evolution of star formation in clusters is the work of 
iPoggianti etail (120061) . These authors used the [O II]A3727 
line as a tracer of star formation, rather than color, and mea- 
sured the fraction of galaxies with [O II] emission (equiva- 
lent width > 3A) as a function of both cluster redshift and 
cluster velocity dispersion. Their sample includes 25 clus- 
ters with z = 0.4-0.8 and another 10 groups in the same red- 
shift range, while they have a large local comparison sample 
at z = 0.04-0.08 from the Sloan Digital Sky Survey. They 
measure the [O II] fraction f[oa] relative to an evolving abso- 
lute magnitude limit Mv,Um that varies from -20.5 at z = 0.8 to 
-20.1 at z = 0.4, while the local limit was M v < -19.8. Their 
main results are that there is substantial evolution in /joii] and 
that there is substantial variation in /[on] with velocity dis- 
persion at a given redshift. Given the velocity dispersion de- 
pendence, a direct comparison of the evolution of f[anj with 
/a is not meaningful for different cluster samples. Instead, 
we have used their upper envelope for /[oii](c) at high red- 
shift and their envelope prescription at low redshift to estimate 
/[Oil] for each of our clusters and then computed the average 
/[Oil] for each of the subsamples shown in Table [8] These re- 
lations predict an increase in /oil] of less than a factor of two 
from the low-redshift to the high-redshift subsamples, or sub- 
stantially less than the factor of eight we observe for the AGN 
fraction. 

These results are interesting, although numerous caveats 
forestall too much interpretation of the relative rates of evo- 
lution. One major concern is that there is likely down- 
sizin g in clusters similar to what is observed in the field 
(e.g.lCowie et al.lll996tlH asinger et al. 2005; Silverm an et al.l 
120081: 1 Yencho et al.l2009l) . that is the relative number of galax- 
ies with star formation or AGN activity above a certain thresh- 
old varies with redshift. The direct implication of this for 
the AGN fraction is that the evolution of the AGN fraction 
over a given redshift range is expected to depend on lumi- 
nosity, just as the rate of evolution of the AGN space den- 
sity is observed to vary in the field as a function of minimum 
luminosity (e.g. lUeda et al.ll2003l) . This is similarly a com- 
plication for interpretation of the evolution of star formation, 
and consequently limits direct comparison of the mere rates 
of evolution of star forma tion and AGN above s ome thresh- 
old. For example, while Poggiant i et al.l (|2006) have simi- 



larly used an evolving galaxy luminosity threshold to charac- 
terize the evolution of the star-forming galaxy fraction, their 
galaxy luminosity threshold is over a magnitude fainter and 
therefore they have measured the evolution in a population 
that includes many more fainter cluster members. However, 
these concerns are not an obvious limitation to comparisons 
that use the same luminosity threshold to separately compare 
either AGN or star formation across different environments, 
particularly when the evolution of the star formation rate and 
AGN luminosity are tied to the same galaxy population. For 
example, if the relative rates of evolution of AGN and star 
formation in < + 1 galaxies were different in the field and 
clusters, this would suggest a limit to the extent of the appar- 
ent coevolution of black holes and galaxies in at least one of 
these environments. 

Another concern about a direct comparison to these mea- 
surements of the evolution of the star forming galaxy popu- 
lation is that [O II] emission is more susceptible to redden- 
ing and metallicity eff ects relative to othe r star formation in- 
dicators, such as H a dKewlev et al J 120041) . Many ISO stud- 
ies (summarized by iMetcalfe et al.ll2005l) found evidence for 
an increase in star formation in clusters at higher redshifts, 
and that the increase appeared to be greater than that pre- 
dicted by UV continuum or visible-wavelength spectroscopic 
diagnostics. Spitzer observations of clusters have also found 
substanti al, often obscured , star formation in high-redshift 
cluste rs jGeach et alj|2006t iMarcillac etail 120071 ; iBai et all 
120071) . iGeach et al.l d2006l) used new Spitzer data for two 
clusters and data for five others from the literature to esti- 
mate the star formation rate normalized by the cluster mass. 
They find evidence for an increase in higher redshift clus- 
ters, but also substantial varia t ion be tween clusters at the 
same redshift. Saintonge et al. (2008) used a larger sample 
of eight clusters with 24/im Spitzer data to study the evolu- 
tion of the fraction of obscured star-forming galaxies from 
Z = 0.02-0.83. They find that the fraction of cluster galaxies 
with star formation rates above 5 M Q yr _1 increases from 3% 
at z = 0.02 to 13% at z = 0.83 and that this is stronger evolu- 
tion than exhibited by color-sele ction, such as the criteria of 
IButcher & Oemleri d 19781 119841) . The star-forming galaxies 
they identify in these clusters are also mostly disjoint from the 
Butcher-Oemler galaxies and consequently when they sum 
the blue and mid-infrared galaxies the fraction of star-forming 
galaxies increases to ~ 23% at high redshift. 

Several of these Spitzer studies overlap clusters that are 
also in our sample and it is interesting to see if there is a 
direct correspondence between the AGN and mid-infrared 
sources detecte d by Spitzer. The massive cluster MS 1054-03 
was studied by Bai et al. (2003> and t heir 24/im sou r ces in - 
clude the two X - ray AG N identified bv lJohnson et al.l (12003b . 
ISaintonge et al.l (120081) have three clusters in common with 
our sample: MS0451.6-0305,MS2053.7-0449,andMS 1054- 
03, although they do not provide information on individual 
sources. While not i n our sample, the study of RX J0152.7- 
1357 (z = 0.831) by IMarcillac et all (f2007) found that the 
two most luminous 24pm sources (of 22 confirmed mem- 
bers) were also X-ray AGN. Similarly, IGeach et all (120091) 
found that one (of 12) of the luminous infrared galaxies 
(Lir > 10 n L Q ) in CL0024+16 (z = 0.4) was obviously an 
AGN based on their i nfrared data alone. At lower redshifts, 
Galla gher et al.l (120081) have also used Spitzer data to identify 
A GN and star form ing galaxies in Hickson Compact Groups. 

ISaintonge et al.l (|2008) explore whether or not the increase 
in the fraction of obscured star formation in high-redshift 
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clusters is related to infall. They speculate that the increase 
in star formation reflects the infall of new members and note 
that most of the MlPS-detected cluster galaxies are not pro- 
jected onto the cluster core (inner 500pc). Over larger scales 
the work of iGallazzi et alJ (120091) explored the obscured star 
formation fraction as a function of environment in the Abell 
901/902 supercluster at z = 0.165. They find more obscured 
star formation at intermediate densities than in the cluster 
cores, similar to the distrib ution of the AGN population stud- 
ied by Gilmour et alJ d2007l) in the same supercluster. If there 
is a substantial increase in the obscured star formation fraction 
in the intermediate densities around clusters, and the star for- 
mation in this environment increases with redshift, then pro- 
jection of some of these structures onto the cluster core may 
contaminate the cluster estimates. 

As discussed in Section [431 AGN in the large-scale envi- 
ronments around massive clusters may also project onto clus- 
ter cores. To better evaluate this possibility, it is useful to 
both directly measure the AGN population immediately out- 
side clusters and measure the AGN po pulation in intermedi- 
ate densities more generally. Just as iPoggianti et alJ (2006) 
found that the fraction of [O II] -emitting galaxies increases in 
lower velocity dispersion environments, the AGN fraction as 
a function of environment is important because the environ- 
mental dependence may provide new information on the pro- 
cesses that drive AGN ev olution. Both the XMM observations 
of the COSMOS fields (ISilverman et al.ll2009allbh and Chan- 
dra observations of the Extended Groth Strip from DEEP2 
(the All-wavel ength Extended Groth stri p International Sur- 
vey, A EGIS; IGeorgakakis et all l2008at iGeorgakakis et alJ 
2008b) have estimated the AGN fraction in groups of galax- 
ies or as a function of lo cal overdensity at high redshifts. 
Georgakak is et alJ (l2008bl) found that X-ray AGN are more 
frequently found in groups than in the field, which they con- 
nect to their observation that the X-ray AGN host galaxies 
are often red, luminous galaxies that tend to reside in denser 
environments, although they also find that this trend may 
reverse for the most powerful AGN. In a narrower redshift 
range from 0.7 < z < 0.9 and for Mb < -20 mag they find 
that the AGN fraction is comparable in groups and the field 
and about 5%. This is approximately a factor of five higher 
than we find in clusters at similar redshifts , although these 
values are not exactly comparable as the IGeorgakakis et alJ 
( 2008b]) AGN include somewhat lower-luminosity sources 
than our sample and the host galaxy magn itude limits are 
somewhat different. ISilverman et alJ d2009al) also investi- 
gate the environment dependence of X-ray AGN hosted by 
galaxies above a fixed stellar mass and find no strong prefer- 
ence between the fie ld and groups exce pt for the most mas- 
sive galaxies, while iJeltema et alJ (120071) find that the frac- 
tion of [O II] -emitting galaxies in intermediate-redshift, X- 
ray-selected groups (0.2 < z < 0.6) is similar to clusters at the 
same redshift. 

The clustering analysis by lCoil et all (120091) on the AEGIS 
data also helps to elucidate the distribution of AGN at high 
redshift as a function of environment, AGN luminosity, and 
host galaxy mass. They find that the X-ray AGN have sim- 
ilar clustering to luminous red galaxies and are more likely 
to reside in groups, while UV-bright QSOs are less strongly 
clustered and more similar to the field blu e galaxy popula- 
tion. T his is also similar to the results from lKauffmann et alJ 
(2004) at low redshifts from SDSS, who find that galaxies at a 
fixed stellar mass that host luminous [O III] emission are twice 
as common in low-density regions as high. Taken together, 



the AEGIS and COSMOS results illustrate that the measured 
AGN fraction depends on both the stellar mass (or luminos- 
ity) of the galaxy population and the star formation rate of the 
host, in addition to the AGN luminosity. This makes a di- 
rect comparison between these two surveys, as well as to our 
work on high-redshift clust ers, somewhat p r oblema tic. The 
X-ray range considered by ISilverman et al.1 d2009al) extends 
over 42 < logLo.5-iokeV < 43.7, or approximately half an or- 
der of magnitude below our X-ray threshold for a ty pical AGN 
SED. The X-ray AGN studied bv lCoil et al] d2009l) extend an 
order of magnitude fainter than our work to a hard band limit 
of Lx,h > 10 42 erg s" 1 . Both of these surveys are therefore 
dominated by intrinsically less luminous objects. The galaxy 
mass and luminosity ranges are similarly not identical. In fu- 
ture work we hope to put all of these high-redshift measure- 
ments on an equal basis for a more direct comparison. 

While none of these results suggest that there are more lu- 
minous AGN in clusters than groups or the fields out to z ~ 1 , 
such a trend may be seen at yet higher redshits. Observa- 
tions of cluster galaxies, particularly massive cluster ellipti- 
cals, suggest that most of their stars formed earlier than fiel d 
galaxies (by 0.4 Gyr; Ivan Dokkum & van derMarell l2007). 
If the central black holes of these galaxies grew contem- 
poraneously, then perhaps by z ~ 2 the AGN fraction will 
be higher in denser environments. Some interesting sup- 
port for this picture comes from Chandra observation s of the 
SSA22 protocluster at z = 3.09 dLehmer et alJl2009t) . They 
find a slightly higher AGN fraction in Lyman Break and 
Lya-emitting galaxies in the protocluster compared to the 
field. While this is just one region, observations of the AGN 
fraction in clusters relative to the field at z ~ 2 and above could 
provide interesting new insights into the coevolution of black 
holes and galaxies. 

7. SUMMARY 

We have conducted an expanded survey to identify lumi- 
nous Lxji > 10 43 erg s -1 AGN in clusters of galaxies from 
z ~ 0.05 to z ~ 1.3. At low redshifts we have presented a new 
X-ray analysis of archival Chandra observations and spec- 
troscopic follow-up of AGN candidates in six new clusters. 
There are no new, luminous AGN in these clusters and there 
are a total of just two luminous AGN in our sample of 17 
clusters with z < 0.4. These measurements further strengthen 
the evidence for a very small luminous AGN fraction in low- 
redshift clusters. An important virtue of the new clusters is 
that the X-ray and spectroscopic coverage extends to the pro- 
jected Z?200 radius and therefore they are better matched to 
observations of high-redshift clusters. At higher redshifts we 
have combined our previous work with literature data on X- 
ray sources, primarily from the ChaMP and SEXSI surveys, 
to compile a total sample of 15 clusters at z > 0.4. In spite 
of somewhat incomplete spectroscopic coverage of the X-ray 
sources in these fields, there are 18 luminous AGN in these 
clusters. 

We parameterize the evolution of the AGN population in 
clusters in terms of the fraction of luminous galaxies that host 
AGN above our luminosity threshold. We have used a vari- 
ety of techniques to estimate the number of luminous galax- 
ies, defined to have Mr < M'^+ 1, in these clusters and cal- 
culated the average cluster AGN fraction in several redshift 
bins. As the low and high-redshift clusters are reasonably 
well matched in terms of cluster velocity dispersion and X- 
ray temperature, the increase in the number of AGN is closely 
related to the increase in the fraction of galaxies more lumi- 
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nous than M% + 1 . Specifically, we find that the AGN fraction 
increases by approximately a factor of eight from z ~ 0.2 to 
Z ~ 1. This corresponds to an increase in the AGN popula- 
tion that scales as (1 +z) 5 3 . If the radio AGN population in 
clusters increases by a comparable amount, radio AGN may 
impact the identification of clusters as a function of redshift in 
current and planned SZ surveys. The substantial evolution in 
the cluster AGN population is also correlated with the evolu- 
tion of the fraction of star-forming galaxies in clusters known 
as the Butcher-Oemler effect. Detailed studies of star forma- 
tion and AGN in individual clusters could better quantify the 
extent that these two phenomena are coupled in clusters or 
perhaps even individual galaxies. We have also estimated the 
evolution of the field AGN fraction to compare it to the clus- 
ter AGN fraction. While the field AGN fraction is higher at 
all redshifts, the present data do not suffice to conclude if the 
rate of evolution is faster or slower in clusters. Future mea- 
surements of the relative evolution of star formation and black 
hole growth in clusters and the field could be an important 
probe of the coevolution of black holes and their host galax- 
ies. 

Measurements of the radial distribution of the cluster AGN 
provide new information on the origin of AGN within clus- 
ters. Unlike we found in previous work at low redshifts, the 
AGN in these high-redshift clusters are not strongly centrally 
concentrated when their distribution is plotted normalized to 
the /?200 radius. This demonstrates that there are substantial 
numbers in the outskirts of clusters and supports the hypoth- 
esis that some cluster AGN are hosted by relatively gas-rich 
galaxies that have recently entered the cluster potential. While 
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